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BENDING MOLES, AND CONVERSION MODE OPERATION 


Wayne Johnson* 

U.S. Army Air Mobility RAD Laboratory 
Moffett Field, California 


SUMMARY 

An analytical model is developed for proprotor aircraft dynamics. 

The rotor model includes coupled flap- lag bending modes, and blade torsion 
degrees of freedom. The rotor aerodynamic model is generally valid for high 
and low inflow, and for axial and nonaxial flight. For the rotor support, 
a cantilever wing is considered; incorporation of a more general support 
with this rotor model will be a straight-forward matter. 


INTRODUCTION 

This report presents the development of an analytical model for tilting 
proprotor aircraft dynamics. The emphasis in this model is on the rotor. 

The rotor support for the present is limited to a cantilever wing, but 
the incorporation of a more general support model with this rotor model 
will be a straight-forward matter. 

The rotor motion is represented by; coupled flap and lag bending 
modes; rigid pitch (control system flexibility) and blade elastic torsion 
deflection; gimbal tilt and rotor speed perturbation degrees of freedom 
(optional). The six components of shaft linear and angular motion are 
Included, and rotor blade pitch oontrol. The rotor aerodynamic model la 
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generally valid for high and low inflow, and for axial and nonaxial flight. 

The effects of compressibility and static stall may be included, but 
reverse flow and unsteady wake aerodynamic interference effects are 
neglected. Three components of aerodynamic gust are included as external 
excitation. The rotor model Includes gimbal undersling, torque offset, 
precone, droop, sweep, and feathering axis offset (for the case with blade 
bending flexibility inboard of the pitch bearing). Center of gravity, 
aerodynamic center, and tension center offsets are included; but the elastic 
axis is assumed to be a straight line, and only offset from the pitch 
axis by the droop and sweep rotations. Kor the equations of motion in 
the nonrotating frame it is assumed the rotor has three or more blades. 

The equations of motion are derived for the rotor degrees of freedom, 
along with the forces and moments acting on the hub. 

This rotor model may be coupled with any support model. The present 
derivation is restricted to a cantilever wing support (fig. l). The wing motion 
is represented by three degrees of freedom; wing vertical bending, wing 
chordwise bending, and wing torsion. Wing aerodynamic forces are included, 
and a wing trailing edge flap in the controls. 

The differential equations of motion for the pro pro tor and 
support system are presented in matrix form, for three cases; axial flow, 
which is a constant coefficient system; nonaxial flow, which is properly 
a periodic coefficient system; and a constant coefficient approximation 
for the nonaxial flow equations, using the mean of the coefficients in the 
nonrotating frame. The axial flow case is applicable to the pro pro tor 
aircraft in airplane mode cruise and in helioopter node hover flight. 

The nonaxial flow case Is applicable to helioopter mode forward flight, 
and to conversion node flight of the pro pro tor aircraft. 

Solutions and results for pro pro tor dynamics from these equations 
are not presented in this report, but are left to a later work. » 
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The body of this report is composed of the following 
sections t 

Bending/Torsion of Highly Twisted Beam 
Equations of Motion for a Rotating Blade 
Aerodynamics 
Rotor Trim 

Blade Bending and Torsion Modes 

Support Equations of Motioni Cantilever Wing 

Equations of Motion 


bending/torsion of highly twisted beam 


This section presents an engineering bean theory model for the 
coupled flap/lag bending and torsion of a rotor blade, with large pitch 
and twist, A high aspect ratio (cf the structural elements) is assumed 
so the beam model Is applicable. The object is to relate the bending 
moments at the section, and the torsion moment, to the blade deflection and 
elastic torsion at that section. The analysis follows the work of 
references 1-3. 

The basic assumptions are i) an elastic axis exists, and the 
undeformed elastic axis is a straight line, and ii) the blade has a 
high aspect ratio (of the structural elements), so engineering beam 
theory applies. Figure 2 shows the geometry of the undeformed blade. 

The span vaiable is r, measured from the center of rotation along the 
straight elastic axis. The section coordinates are x and z, the principle 
axes of the section, with origin at the elastic axis. So by definition 

C ** <JA =0 

Really the integral is over the tension carring elements, i.e. modulus 
weighted, ^xzE dA ■ 0| so x and z are modulus principle axes . 

This remark holds for all the section integrals in this section. The 
tension center (modulus weighted centroid) is x^, aft of the clastic 
axis, and on the x axis, i.e. , 

© 

Again, these are modulus weighted integrals. If E is uniform over the 
section, then x^ is the area centroid) and if the section mass distribution 
is the same as the E distribution, then x^ equals the section center of 
gravity location. 

The angle of the major principle axis (the x axis) with respect 
to the hub plane is 6 . The existence of the elastic axis means that 
elastic twist about the EA occurs without bending) vnt may, and shall, include 
the elastic torsion deflection in ^ . The blade feathering axis (FA) is at 


r FA' the blade pitch is described by root pitch 0° (rigid pitch about 
the FA, including that due to elastic distortion of the control system), 
built in twist 0^ , and elastic torsion about the EA dt * 

9 — 3-v-w -*■ ©e- 

® W} = root pitchi command collective and 

cyclic and control system flexibility} 
rigid pitch about the FA; equals Q at 

r FA 

* built in twist, ©-**<«>* 
m elastic torsion, d«, > ***^ ■ r 0 

There is stress in the blade due to only. It is assumed that Qu is 
small, but 0 and 9-tu are allowed to be large* 

The unit vectors in the hub plane (HP) axis system (rotating) 
are tf® y "ja • The unit vectors for the principle axes of the 

section (x,r,z) are } these are for no bending, but including 

the elastic torsion in the pitch angle © . So the principle unit 

vectors are rotated by 0 from the HP: 

*r l?6«a»sd — Teg 

^ =r -t- 


Description of the bending 

Now the engineering beam theory assumption is introduced} ill) plane 

sections perpendicular to the EA remain so after the bending of the blade. 

Figure 3 shows the geometry of the deformed section. The deformation of 

the blade is described by 

i) deflection of the BA: x-, r , 

0 0 0 

ii) rotation of the section! 
lii) twist about the BA, implicit in ,*£ 

The quantities x q , r Q , s o , ^ are assumed to be small. 
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The unit vectors of the unbent cross section are * The 

unit vectors in the deformed cross section are 1 3^*6 

sure the prlciple axes of the section, and is tangent to the deformed 

EA. It follows then that 

3 t if 

-&.S ~ t - ^ 


Now by definition s {)r/&s where r =■ 4*<«* + *«y5 + 

and s » arc length along the deformed EA. Hence to first order 

-?** * -f +■ 

j -H t%l + (*i -Vo*' 

It follows the rotation of the section is 

— 4>* = >4 + !•©* 

“ !• ” *• ® 

or = c 

The undeflected position of the blade element is 
r— • r^+xT?+ 
and the deflected position 

fc+ivyj + taS + + +"6%% + *H*i 

, r J r 0 -}+ ~*4 >*V? *h*T? + it 



He shall neglect r Q for now. The «W.k analysis is simplified then since 
to first order s - r» r Q just gives uniform strain over the section, so it 
nay be simply added back later. 
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Analysis of strain 

The metric of the undeform erf b]3 da — no bending, and no torsion 
so — is 

r = xV + r*3> + fclc 

|? = -*■ "5 

O'** 

3 rr s d"7-'dr - x 4 ’ ®‘ rw * -> 

The metric of the deformed blade, with bending and torsion, is 

(2 - (*■»■*# VT +1 r +*• 4> % — +■ ^ 

d ' 

o J 

Q rr a ^ “ O •* X <Vl - ■+ -td"*****#^ V’ 

+ - e v < x x 

Then the axial component of the stain tensor i 

5 rr = t «*"- - 3 rr>) 

-*• (*» — ® , (»+X»^V’ — 6w > l 3 

The linear strain, for small x Q , z q , 0^ , • ^ , is 

Krr » 6cr *♦* ©i* 

+ e£, C*4 -*%: +el 
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The strain due to the blade tension, , Is a constant such 

that the tension is given by 

-T sr { £*rr S > A * S E£>A 

Substituting for €rr and using dA ■ 0, ^x dA = x c A, and 

5 + ~ 3^ a li?A 

where kp is the (modulus t-. ghted) radius of gyration about the EA, obtain 
for e T i 

* 4*;j. *0*C “ ©iv» 

-4- & St, -I" Co 

In this expression, the strain due to the blade extension r has been anded, 

o 

It follows the .'.train may be written, wl Lh *v = T/EA i 

<£rr - €.t ■+■ C * — ^ “♦* ©tu> 

•+ dni d c +l l - Vf ) 


Section moments 

To find the moments on the section, the second engineering beam 
theory assumption is introduced: iv) all stresses except *cr are 

negligible. The axial stress is given by • The direction 

of *rr iB A a2/|a£i 

The moment on the deformed crocs section (figure 4) is 

To find M, integrate the moment about the EA due to the elemental force 
v m tr OK on the crocs section: 

as Cut?** + * *?C 

= £-*■£** + -+ ©A 

Integrating over the blade section, there follows the total moments due 
to beading and elastic torsion: 
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Kx = - S ^ ^ 


= 


S 




Tr 


9k 


(*c 


= csadt -»- 5 (V + % 1 *> &-^, 

/*4^hO*S. 


tw N rr 


S>A 


To !•! has been added the tore ion moment GJ ©» , due to snear stresses 

r 

produced by elastic torsion. There moments are about the &•. . For bending 
It is more convenient to work with moments about the tension center x„: 


A\x = — $ 


Substituting for and integrating, the momenta arei 

Mk — £ ( #>*-*• 4 > i ^ — © y 
« EJTu* ( 4 >* — 4 **^ + e> y 

Aa (. s ^3"7 4 T + (3 E c?e, + dtw ^ 

+ e , t£ 3 : xt (>£-’ a * 4 »^ + 0 * 0 ] 

where 

XT* . = 5 *c"> X 0 A 

» Vf A = ^(x t 4*i t ^QA 

.x.^ st ^ (x-Wi ( x 1 -*-* 1 ') &A 

*X%* » $ \(* x + % ty )bk 

^px = 5 (X* + 
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The Integrals are all over the tension carrying elements of course, i.e, 
modulus weighted. The tension T acts at the tension center x„s hence the 
bending moments about the EA are given *'om those about x^ by 

/ Ai OeA - 

The bending/ tors ion coupling is due to EI^p and EI^; ^ or a symmetrical 
section EI^p - 0, 


Vector formulation 

Define the bending moment vector 

and the flap/lag deflection ^ 

( *t is not quite the moment on the section, because and are 
really the *£«t and "W«( components of the moment). The derivatives 
of w are . 

Cfc»T? — *♦* 

Then the result for the bending and torsion moments may be written i 


Mg 1 =• ('ex %% *dT? + • (*.*£ —x.iY" 

4- ©£ c ex«tt -ex«* t } 

M r = (>*3 * kfi'-v + ©i« k|-r- 

+ ©£, c cr^-ex t ^^ . Ct.**? 

Ibis is the result sought in this section* 


- 10 - 



1 


Writing the El dyadic as and ^he 

coupling as £Xj»* , this result becomes 

•J tah J ff y / ^ 

— E^- w ■+• ©^» 0 e Elf 

/v^ r = ^OTT + Vf t + 0twj E X. f>f~^ d*, ■+■ fcf fStvi 

♦ — — * J«y 

-+- ©-pj t -Cf • w 

This form is an obvious extension of the engineering beam theory result 
for uncoupled bending and torsion ( ■ 0 case). The vector formulation 

of the result is a major simplification. The vector form allows an easy 
transformation from one axis system to another. In fact, the vector form 
is independent of the axis system used (the base of the vectors), which 
is the source of the simplification. Working with the vector form simplifies 
the analysis to follow j the base of the vectors (for example, either the 
hub plane system, , or the principle axis system, x? ) 

will be considered only when come to evaluate the coefficients of the 
equations of motion, never in the derivation. 

This is a linearized result. So the ^ appearing in El 
and in”w are based on the trim pitch angle 0*©*-*^*. The perturbation 
of due to gives second order moments, which have already 

been neglected in the derivation. The net torsion modulus is 

* <»-x + 0 ^ ffXf* 

where T • • centrifugal tension in the blade. For 

the elastic torsion stiffness characteristic of rotor blades, the GJ tern 

usually dominates. The k« T tern is only important for very soft (torslonally) 

* ^ 

blades, near the root. The Orw tern is only important for very 

soft, high twist blades. 
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SGp’ATit NS IT-’ ::CTHN i'l A ROTATING BI.AjE 

rhir section derives the -tions r ( t.*n i for a helicopter 
rotor blade, .'he blade ~.otion considered ln< i c coupled flnp/la g 
Pending (including the rigid modes if the blade is articulated), 
rigid pitch, elastic torsion, gimbal pitch and roll, and rotation speed 
perturbation. The analysis includes the effects of precone, droop, and 
sweep; feathering axis offset; and rque offset and ginbal undersling. 

The effects of shaft motion, and the hub forces and moments are Included, 
so this analysis may be combined with the equations of motion for a body 
or support to give the complete aeroelactic model for the system. 

Numerous approximations are made in the course of the analysis, in order 
to obtain a tractable set of equations. 

Rotor i( Confi£i2^ 

Consider an N-b laded rotor, rotating at speed S~2. (figure 5 ). 

The rath rotor blade is at 

H',*, * 

where » 2*it/N and *r» Sit is a nondimenaional time variable. 

The S system (*i* J *A ) is a nonrotating, hub plane coordinate 

system; it is an inertial frame. The B system ( **?§ ^"€4 ) is a 

coordinate frame rotating with the nth blade. The acceleration, angular 
velocity, and angular acceleration of the hub, and the forces and moments 
exerted by the rotor on the hub are defined In the nonrotating HP frame -- 
the S system. The rotor blade equations of motion are derived in the 
rotating frame — the B system. Figure 6(a) shows the definition of the 
rotor shaft motion, linear and angular displacement in an inertial freme. 
Figure 6(b) shows ths definition of forces and moments on the hub, in the 
nonrotating frame. 

Blade root geometry 

Figure v shows the blade root geometry considered (wvx)i*tortt,d). 
The origin of the B system is the location of the glmbali if there is no 
gimbal, this is just ths point where the shaft notion and hub forces are 
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defined. Tne gimbal is at the center cf the B and 5 frames. The huh 
of the rotor is z pA below the gimbal (gimbal undersling). The torque 
offset Xp A is positive in the direction. The azimuth «k* is 

measured to the feathering axis line (its projection in the HP), so the 
FA is parallel to the dk axis, and offset from the center of 
rotation. The precone angle Sp A* gives the orientation of the FA with 

respect to the hub plane; <&*, is positive upward, and is assumed to 
be a small angle* The FA is offset from the center of rotation by r v ; 

- n 

the FA is located at r = r pA along the blade. Th^ rigid pitch rotation 
of the blade about the feathering axis occurs at r pA * The droop angle 
£>PK Z and sweep angle <$9*3 occur at r pA , just outboard of the 
feather bearing; and <^3 give the orientation of the SA of 

the blade with respect to the FA. Note that these angles are measured 
in the HP frame; £pa % is positive downward, and is positive 

aft. Both and £9*3 are assumed to be small angles. 

From the gimbal to the blade root is the hub, underslung by z p ^ 
and torque offset by Xp A « From the root to the FA is a shank, of length 
r FA' undistorted is . straight line an an angle <5pa, to the hub 

plane (small precone). The blade outboard of the FA at r pA , undistorted, 
is a straight elastic axis, with small droop and sweep ( and ) 

with: respect to the FA direction. 

From the gimbal to the root is a rigid hub. The shank (inboard 
of the FA at r pA ) and the blade (outboard of r pA ) are flexible in bending. 

The shank is assumed to be rigid in torsion, i.e. the effect of torsion of 
the hub inboard of the feathering axis is neglected. The 'lade outboard 
of the FA is flexible in torsion as well as bending. There is rigid pitch 
rotation of the blade about the FA, which takes place at r pA , about the local 
direction of the FA at r pA , including the bending of the shank. Inclusion 
of the bending flexibility of the blade inboard of the feathering axis 
means the the general rotor configuration is considered; the articulated 
rotor with the FA inboard or outboard of the hinges, or the cantilever blade 
with or without flexibility inboard of the FA. The special case of a rigid 
shrink can be considered as well of course. 
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Geometry of the blade 

Figure P shows the undeformed geometry of the blade. It is assume' 5 
that i) an elastic axis exists, and the undeformed EA is a straight line; 
and ii) the blade has a high aspect ratio, so engineering bean theory ana 
lifting line theory are applicable. The following notation is used: 


FA 

SA 

CG 

AC 

TC 


feathering axis 
elastic axis 

locus of section center of gravity 
locus of section aerodynamic center 
locus of section tension canter 


The distances y T , x, , and x„ are positive aft, measured from the EA; they 
are in general a function of r. The corresponding z displacements are 
neglected, i.e. taken as zero. 


The system is the EA/principle axis system of the section. 

The subscript "o’* is for the undeformed frame, i.e. no elastic torsion In 
3 , or ginbal or rotor speed degrees of freedom. The subscript will 

be dropped when it is obvious what is meant. The direction of the undefermed 
EA is ; are the directions of the local principle axes of 

the section, undeformed (no bending or torsion). 


The span variable is r, measured from the center of rotation to 
the tip. This variable is dimensionless, r * 0 at the root to r * 1 at 
the tip. The section coordinates x and z are mass principle axes, with 
origin at the EA, It is assumed that the direction of the mass principle 
is the same as the modulus principle axes (used in engineering beam theory 
for the structural moments). The CG is at z - 0, x - Xj. Usually Xj and 
x Q should be close. By definition then 

C s vw section mass 




CG location 


axes 


and 
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section polar moment 
of inertia, about EA 


1 


The blade pitch ancle is 0 ; here undistorted, denoted by the 
subscript "m". The ancle 0 is measured from the HP to the section 
principlejkxis . It is then the angle of rotation of from the 

HP axes.' The undeformed pitch angle is the collective plus the built 
m twist © = 


where 


* collective pitch 
dru = twist 


Define 6^ as the pitch at so 0. The root pitch is 

then ©° = dt»\t . Inboard of do not have the rotation of 

the blade, but there can be pitch of the local prinicple axes with respect 
to the HP, which is included in d-*u» for r < r p A » Note is 

not necessarily zero, hence there is a jump in 0 at r^ of magnitude 

© cr >£ - O C<>* ^ - €>^Cr p “ ^ 

The trim pitch angle is then 

e = e«.= \ 


It is assumed that 8 m is steady, constant in time, so independent 
of . Cyclic variations in 0 , as may be required to trim the rotor, 

are included in the perturbation to the pitch angle. We shall allow 
the trim pitch angle to be large, hence ^ and 0-nj may be largj 
angles. 


The physical sweep and droop angles are defined with respect to 
the blade outboard of the FA, i.e. rotated by 0° about the FA. Let 
and <5**3 be with respect to the principle axes at the 

root (at the FA, r * r^) j these angles will be equivalent to and 

when there is zero root pitch. Hence the droop and sweep angles are 


! ’ 
i J 
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-*-£pa x ^ /v ' e> ° •* o^a© 

The angles <Jp^ t and «&FAj are fixed geometrical constants. It follows 

then that 4fA t and ifA, vary with the root pitch 6° • This must 
be accounted lor when there are perturbations to O due to the rigid 
pitch motion of the blade. In addition, the droop and sweep only affect 
the blade outboard of the FA, i.e. for r> r fA . This may be accounted for 
by including with <5*^ and <SpAj the factor U(r-r pA ), where 

r i r>o 

lA ( <”> — 2 o r 4 o 

rfe shall follow the convention of assuming the factor M is present whenever 
writing or . 

From the B system (mth blade, rotating HP axes) to the o system 
(undistorted EA/XS axes) there is 

1) rotation ip*,- about (small precone and droop) 

2) rotation about kg (small sweep) 

3) then rotation about ICA (large pitch angle) 

Hence 

■£«,* ^ SM 8*. _ £pAj W49- ^ 

“fo* + bSrk^l + 

where and are based on ©£ m ©u»u , and are absent for 

r <• r_, . . We shall drop the subscripts "o" and "m", denoting the trim and 
FA 

undistorted geometry, when it is obvious what is meant. 
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Motion 


The rotor blade motion (degrees of freedom of the rotor) is 
described by: 

1) gimbal motion (optional); pitch and roll of the 
rotor disk. 

2 ) rotor speed perturbation. 

3) Then elastic torsion about the EA, and rigid pitch 
about the FA. 

4) Followed by bending deflection of the FA, including 
rigid flap and lag motion if the blade is articulated. 

Gimbal motlen/rotor speed perturbation 

Figure 9(a) shows the gimbal motion and rotor speed perturbation 
in the nonrotating frame. The gimbal degrees of freedom are and frs ' 
rotation of the rotor disk, in the nonrotating frame (S system), with the 
same convention as (3, c and tip path plane tilt. The rotor rotational 
speed perturbation is . The degree of freedom H> s is a rotation 

about the shaft axis ; so the azimuth angle of the mth blade is really 

•K* + 

Figure 9(b) shows the gimbal motion in the rotating frame. 

The iegrees of freedom are ^ and , given by 

(5<jr * (^<*G t** 4- ^Qrjt K.H’*. 

■+*■ (3<>f co sVm 

The main effects are due to , the flapriae rotation about the Tfg axis; 

, the rotation about , only introduces a translation of the 
hub due to Sp A and x^. The blade pitch © is defined with respect to 
the hub plane, so only the blade inboard of the FA sees the pitch rotation 
due to , and that effect will be neglected. 

Blade motion 

Figure 3 shows the geometry of the deformed blade. The blade 
deformation is described byt 
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l) twist about the EA: & 

?) deflection of the EA: x , z 

o o 

3) notation of the section: 


The pitch angle & , including perturbations, is implicit in 

the systemi tit are the principle axes of the blade 

with no bending, but now with the blade elastic torsion and rigid pitch 
motion in © .The XS system ( } ) are the principle 

axes and Ea of the deformed blade, including torsion and bending. The 
tangent to the deformed SA is ; the rotation of the cross section 

from is given by and 5 


The blade position, relative the root, is then: 

^ — f f + r < o'y-f + -v vc* 4 + %Vl* s 

~ + (wc+%«iD -f 

^ivy V 


% \ 
\ A 


°A 

\v 

V 


> 5 .\ 


"i 'A 


We will neglect the perturbation of the radial position, f.-f 


Blade pitch 

The angle is the angle of the major principle axis of the 

section (the x axis, chordwise), measured from the hub plane. The blade 
pitch is composed of: 

1) e*(H^ “ root pitch, the pitch of the blade at the FA 

at r - rf. 1 due to commanded collective and control, 
control system flexibility, and mechanical feedback. 

2) ■ built-in twist; 0. 

3) ” elastic torsion about the EA; aero at the FA, 

r>* ^ “ 0, only produces torsion shear 
stress in the blade. 
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For the shank, r^r^, elastic torsion is neglected, and it does not 
see the root pitch . Then is used for the pitch of the principle 

axes with respect to the hub plane in the shank. There is no perturbation 
to © inboard of the FA, the pitch and torsion degrees of freedom are 
only for outboard of the FA. Since probably is not zero, there 

is a jump in © at the FA. So the blade pitch is 




The commanded root pitch angle is 

© — 1 + ©TAA 

where 


© = 



dcau = collective pitch angle; the trim value, which 

may be large but is assumed to be steady in time. 

0c** * control input; time dependent, but assumed to be 
a small angle; includes cyclic to trim the rotor; 
and for dynamics analyses this is the control 
variable. 

The blade root pitch commanded by the control system is © c ; ©* is 

the actual blade root pitch. The difference t © fc } is the rigid 
pitch motion due to control system flexibility or mechanical coupling in 
the control system (i.e. effects). Hence we may write the blade 

pitch ast 


( C&um\s + Br*') + d©°— -f ©caw -V © t 


T><>* 


rar, 






Now the pitch © may be separated into trim and perturbation 


C ©* ♦ © 

L 


contributions : 
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where the trim terms are (as above) 


= 


(" d t/>\\ 

0-W* 


■+• €*tv> 


and the per t>irbi t ions 


© 


( q °~ e c> ) -*■ 9*. 

©° « C9‘- © C> > 4* ©, 
O 


The trim value of the pitch is , composed of dc*\t and d-n* ? 

it is a large, steady angle* The perturbation of the pitch angle 
is & , composed of the blade motion . and ©*, ; 

all are small angles, so S is small* For the rigid pitch degree of 
freedom we shall use p o , defined as 

= e° = c ©° - e* ) 9wv. 

and for the elastic pitch ©t* an expansion as a series in the normal 
modes (described in the sections to follow). Note that p Q is the total 
rigid pitch perturbation, including the control 
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Coord^r^e^ranes^an 

S System : nonrotating, hub plane frame 

rotation V*.- ‘ 

B system: rotating, mth blade, hub plane 


H system: hub frame 


FA system: blade FA (EA for r <- r^) 


EA system: EA outboard of PA 


about ^ 

A, 


X 

abo tt 

% ' 

©<v about 


H'j about 


about ~£ H 



— abou* ^FA < 

- SFAj about V-Pa « 


A aDout 
— about 


blade system: principle axes, including torsion 

XS system: principle axes, torsion and bending 

B system 


about 

J>, about Tc 


ft** "^S 
c^aHU 

t* 


- c-oaV^. 


Blade system 

From the B system to the blade system, there is first rotation 
about and rotation %- about i then 

rotation 0 about C*«A . Hence i 


-FI- 
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~t = -£ a — 

£ C^Cr if a, - <fa x ^ S^S ■+• iPAi^ { * J ® 3 

”£ a s^~© "^a ■+■ CjoA 0 Ita 

L- <(3<i + i?*, — SPA^***© + C'+'j- Sf*j^Sn~&3 
M 3 *£^ ~ ~ C ^5 — ■** ((i& — 

For r < rp A * the £* k z and $«k s terms drop; in particulars 

= T^A = "*4 + C (*« **■ 

XS system 

•f us - "t <J>* ^ 

"txs * *fc - ♦» j 

Undisturbed blade system 

The undisturbed blade system is ~€y£)^ without ( 3 $ , , 

or the pitch perturbations in 6 (and ^ <Wj based on 0^39^ ); 

hence « 

**?• * - W- 9 -M -► "tftL <*•*•“ -Jsa, 1 

l?*a *~©~ 1» teee*.*^a + -ft (j^*, - $f Aj s~©~l 

“if* ■ + iwj^t + ( ipx, — Vfa0*^4 

Now since the blade motion 8 , <*• , and H*s is small, it is possible 
to expand the blade system in terms of the undisturbed frame t 

-t me [_ ($+ - - 4 *( , N 4 4 '© # 

”£ * * 4 - © -f :?• 4 > 0 % 4 -S*i»A t 3 &%-«e 3 

There follows them 
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f 


T 


*0^ =■ 4 %®v.® ' -1 © ( ^® 

which is an expansion of the bending/ torsion deflection of the blade 
In terms of the undisturbed axis system. 

Rlade position, velocity, and acceleration 


Position 

The distance from the gimbal to a point on the blade section is 

"<r— - ~ *$*■£»* ■**$* ^$*— '^«A^4f'^4Xo£4*«£4*'£ li +*iL*j 



which may be written 

*** — 

4 3*» (. ~ H > *'> 

4 *^s C — 4 **A ^FX* .) 

4 4 Cx.-e 4 !•*& } 4 ( %-t 4 ^ 

* (-*FA -If*©* -«*■#» 4 Cr-r'i^&n 3 } 

4 3* c r 4 «V* - x* * *Vj ) 

4^ (— %** 4 *** 0* 4 * if a.^ - ( r ~ CV*>y SfA*^ 

4 <*•*£ 4 4 C**£ 4 ) 
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Velocity 

The velocity of a point on the blade, relative to the rotating 
frame (the B system) 1st 

“ ( Sr ^ )ft “ ( - =*Fa ©<* _ r HO 

“** "3*0 ( *FA — *^A H's ") 

( *pa ©& ■+* r ") 

-t-(.r..r PA ^ ©• C —X n 6*k z -T* &? 

■+■ < (Xo + x'iT? + (*• + 3* 

where 

(c**+*vt 1**** 'it y s ( *,-t 0 + i.'&y 

+ © U^• + ^')-L o -(».+’0'£©') 

Acceleration 

l^e acceleration of a point on the blade, relative to the rotating 
frame, and neglecting the squares of velocities, isi 

a r ~ Vto Vf = “to ( — ^?a ~ r H*t ) 

( **A ^ - **A »iO 
*»-^6 <*?a ©«» ^ 

C r — < 9 A ^ e° C — ^^Aj^ 

+ <(*•+»>*£ -+ (*. y 

where 

COo+*')‘i + s C*e*t ••*•!•*£• V* 

-f* 0 ^ 

Acceleration of the blade 

The acceleration of the blade Is required with respect to an 
Inertia frame, l.e^ln the S system. The B system rotates at a constant 
angular velocity with respect to the S frame. The shaft j 

motion la composed of linear and angular velocity and acceleration of 

I 

I ■ I 


. \ 


\ 



the origin of' the S frame (the gimbal jr>int at the hub center of rotation). 

The acceleration, angular velocity, and angular acceleration of the o 

system, with respect to the nonrotating, inertial frame, are: 

» .. •• -5 

do * x* "^5 ■+ 3* + w Ws 

U^o * -* ^4 

^ "** 

_j — b JLi 

It Is assumed that a Q , u>® , and u. are all small quantities. 


Given above is the motion of the blade in the B frame, the acceleration 


and velocity of the blade a and v 

r r 


Now we shall derive the acceleration 

of a blade point in inertial spac* ( "5* ), in terms of the motion of the shaft, 
the rotation of the rotor, and the blade motion in the B frame. From the 
result for the acceleration in a rotating coordinate frame, there follows: 

o = a, -+■ -f- -t- u»« * c u»#' r 

where a and v are the acceleration and velocity of a point in the 
r,s r»s _j 

S frame. The B system rotates at angular velocity it • St Kg with 
respect to the S frame. Hence with _T2. constant and no angular acceleration 
or acceleration of B with respect to 5, there follows: 


o c^s 


or -4* 2 SZ. a vr -+■ St M C 5Z x ) 


^ O i Q 

Vf^s ~ Vr “t" SZ X 


where a^ and v r are the acceleration and velocity in the B frame. Thus: 

•J - «• -4 2 x w r -+ & * ( Stir} 


-f 2 w»«x 


"^ f 2 G*« x ( JL x ^ ") + vj* x( 4- x** 


To first order in the velocity and angular velocity, this becomes: 

-a -> ^ 


— a a. 
a s o* *r 


ff ■+• 2 Jl » y r a <( Slur ^ 


. j - -> a\ 

-+■ 2 ***• * C sz, xc ; 
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The six terms are respectively: the acceleration of the or.' gin; the relative 

acceleration in the rotating frame: the relative coriolis acceleration; 
the centrifugal acceleration; the coriolis acceleration due to the angular 
velocity of the origin; and the angular acceleration of the origin. In 
dyadic operator form, and with SX. * 51 V-a , this result is 

-+ 2 St C Eg O v r — SZ? -t 3*^4^ r" 

2SZ C 

To obtain the forces and moments and equations, of notion, the 
acceleration is multiplied by the density of the blade point ( dndr ) and 
integrated over the volume of the blade, to produce the total acceleration 
of the blade. 


Equations of Motion and F orces 


The equations of motion for elastic bending, torsion, and rijig pitch 
of the blade are obtained from equilibrium of inertial, aerodynamic, and 
elastic moments on the portion of the blade outboard of r: 

-W*i 


where Mg - structural moment on deformed cross section, on the 

inboard face; so - M_ is the external force on the 
outboard face, 

“ total aerodynamic force on blade surface outboard of r. 
Mj * total acceleration of the blade outboard of r. 


Kg is a general elastic constraint: froa engineering beam theory for 

bar ding and torsion: from control system flexibility for rigid pitch: 
hub spring for gimbal notion: or it is the fores or moment on the hub 
due to the rotor (so - Mg is the fores on the rotor). Hj is the angular 
acceleration of the blade outboard of r, about the point ? 0 (r)> 


Av 


\ 

IrO 
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tor bending, engineering beam theory give 6 

1*2 ^ m. /*<"? -+ = Ct't'H ■+ V.V.*On*£ 

— ^ i 

So this operator is applied to and also. For bending the moments 
about the tension center x p are required. Then the desired PDE for bending 
is obtained from fA * 

For elastic torsion, engineering beam theory gives - * 3 ** -Mg 
So this same operator is applied to and M^. For torsion require moments 
about the section EA (x - 0) at r; also, elastic torsion involves only 
r> r PA* The desired PEE for torsion is then obtained from ~ r 

The equation of motion for rigid pitch degree of freedom p Q = © 

is obtained from equilibrium of moments about the FAi 

v -J> * 

- -3SS ) * M 0>A> 

A 

where M Is the moment about the FA (x * 0) at r ■ r^. The elastic restraint 
from the control system flexibility gives the restoring moment about the 
FA, completing the desired equation of mow. ■>,. 

The equations of motion for the gimbal degrees of freedom and (£$* 

are obtained from equilibrium of moments about the gimbal: 

Mh — 1?^ ‘ M 

i*) - ~i % • ^ 

where M is the total moments (from all N blades) about the gimbal 
point, in the nonrotating frame. 

The equation of motion for the speed perturbation degree of 

. -j 

freedom 7s I s obtained from equilibrium of torque moments Q * ** * t* • M 

a* 

where again M is the total moment about the gimbal point. 

The total rotor force and moment on the hub (at the gimbal point) 
are obtained from a sum over the N blades of F and M , the force 
and moment due to the mth blade t 


- 27 - 


5 = - 


-C. 

vm a t 






-J (W') — > /'vm > 

Since - F and - M 


Z. M 

V"— » ^ 


(v-> 


are the forces on the b 1 ad e , there follows from 
force and moment equilibrium of the entire blade t 


— e 


l V* > 




-> o 
- Fx 

N'a - 


The hub force and moment are required in the nonrotating hub plane frame 
(the S system); the components are defined as; 

f - H-e* 

Note M produces the gimbal and rotor speed perturbation motion, if those 
degrees of freedom are used, but it is also transmitted throught the gimbal 
to the helicopter body or support. 


Aerodynamics 

The aerodynamic forces and moments on the blade are obtained 
from the integral over the span of the aerodynamic forces and pitch moments 
on the blade section. The aerodynamic forces and moment on the section 
are i 

P in hub plane, positive in drag direction, 

(Z* direction), at the EA 

F normal to the hub plane, positive up (kg 

direction), at the EA 

F radial, positive outward ( S’g direction), 

at the EA 

H & moment about the BA, positive nose up 
The forces on the section are F, F, and F j these are the component 

X Z T 

of the aerodynamic lift and drag forces in the huh plane axis system 

(the B frame). F y is here Just the radial drag force; the radial components 
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1 


°f F x and * z due to tilt of the blade when it is bent are included 
exjj.icitly in the results below. 

The aerodynamic force on the section, at the deforced EA, including 
he effect of the rotation of the section due to bending, is thus: 

*+■ 


=■ ■+ 


- ^4r ( P* + &pa, -^fAj )) 


and the aerodynamic momenti 




Equations of Motion and Hub Forces/Moments 


Bend lag 


The equation of notion comes fron 


S- 

dr* 




*7* 


where M is the moment about the tension center (x = x n ) at r, and 




M — * ^ ^ *.o~t + ^ ^ W 


Inertia: Considering first r^r^^, the moment is 


—3 

^x - 


- "W ) »« 


= 

s:sp- 

O ^ *+• -t* (%# + “1 

)C"cf &A** SL^ 

_D 

__ 

dr “ 


i>o-t -V -V *C> y } * \ r ^ 




( y "u + “ Xoi ^ X « 


a? 

- 

C « 8 a»» 

^ J 

- £ C*o“i *«■ fefc + XcO* ^ 

« SUa»Qo *3 


— 

£ ^CxT? < 8 *^* 



(^•+* 3 t — (u. + xVi 

fro't —%<ri - *«?£) \r . 

| *c? {&*** 
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So 


WV 


b % — 1 ^ 

5r’- Mx = + 

— ^ X ^ GK £>~ 

-4- C *~y 

’*’ Lc^-e -x. - ^ -*(£'? ^ j v 

f “I v ,** 

*“ ° XL ~ J * J 

shall neglect the last tern in this resalt, U*4 + T -£ V' 
as order (c/R) smaller than the first term. Including the * 

“ hlCn ^ ““ “ — - -op and sweep, £ 2^^' 

J~fT- — S*. ^iUv 

-*■ f $(**-*$ + *c.Tt )-J -« 

-t [(«.<-*.■*- *.#y ^ s ■}■*?*-««]* 

— i<e- f«,') •'’SfA l ‘L«-#-4(» 3 'ti) A**S^ 

Kh,r ’ *° ls the delto Motion, i.e. „ ^3, at r|to 

a) shaft motion i with r»rj> a ^ve 

* = *• _ ?*£?. 

= 2SCr- -j> 4 - w. - r ^gx'wp 

So 

» «-*•»« = - 

H- 2-fc.r^ < Tfa'^a') ut 
4-^ CTiati + u c 
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b) 


relative acceleration i 

— ~ T* % Su- 

^ <"* ^ 

r -> 

— k« C — ^pa 



-+-T& C V-FA. &a + (i a 
-f- C^oTf-x^V' 

— 0 (^ <* .» + *X ^ "C + 

_e° 


c) centrifugal acceleration: 


“i = - sx l + i r 

— St 1 Tc* (_-*fa — r PA -»-Ta • /J 

-+• .SI*' r / (i<* 4 -i, -Sj^ 

• <* = — Jx z r 

30 r I 

= -sz i I 

I + [(nc>'Odi^T - 

I -ifr-rp^ ©• (Jjtg- J 2 TcO sVta ^°5 

-*©° -[>*,?* -r FA ^ti 3 

L +'S* ***?*©& — 

_si l r KwoSf-i’ - 1 cx.sv C y ap” 

1 Cixti + 4 jT 2 ») l^ A 

I -Kj^Sj — x x «*»e] 

L w>r 
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|XVf 


d ) coriolis acceleration! ci — ? 

5o ISl J -r4j - tVCtet-Xo'* V } 

* 2 S? | 3 * 6 -^V 


V. ^ 

■+■ <E + [— r ^i ~K| (i#i — VO 


For ^.•v r , it is here necessary to include the effect of the change 


in th^-adL-a position of the blade due to bending! 

c wc + tf Y - ( % - V> fa 




*£[ 


+ A| — 


f *4 


Vr — — So C ” Va "^’^ • C 

- te.*e ->.tv . c a,- ii - s 3 ^ 

— (*<• To * ( *•“£ -Xo^ — -(t- c ?y)&2'3 

— ^ -«-*£&- (*#-£- r - (r- r pA ^^3 3 


iV«* 


2 S 2 . 


So ^•Tf-Xo^ . (%.T?-V£-*x'£V' <&^ 

— ^4 C-?«?a 4 CE^l.Tf-^.f-Xtf )^rA,-(r-r^S i 'j 

— ^ H* [_*PA + 7 t|| • ( - y*T? -*x^ -V r fA S3 3 

~ Oc-^t- (**s V > 3 * 

-[(ve-»^-» c t^ 5 f C^j +t* )~o^] T 

■+ S< c ->a^ W?» ■*■ Sj'tg’S Sr PA fs^‘ -*•£ Y ) 

U-<A 1 - 4 ,')^x Crt.t-».tv*e«f, T*] 


-? 3 ~ 


\ 



elastic : 


. . , . -|*v 

L?* a f J 

* Cl „ 

“*" — EX^T?)©** 3 

aerodynamic > The moment about the tension center (x = x^) at r d’ie to 
the blade loading at the EA at i 

M/K — \j^oe " ^ W* t ® ^ * ^CXA-U» ^ 

* S) (*tu 

s» 

v ? 0(1") —3 _ -S pn O 

rJZ* * -f «?«*-> - Fj<-» -v; t* 
art 


Elastic torsion 

The equation of motion is obtained from 

— £ Mr E - ^ r * - ~ 

where ft is the moment about the EA ax r, and 

d ^ _ d_ 4 4 ..ji s -f . -4- f C w*d -*■ m "3 

*r - J ar L 

inertiai 

*«- Sr!)w~. C *V " 

0 C C (A-^ ^ "Iw^SLy-fiO 

* * r ) L ^ J 3 

So ££ — — (-J -► <»rf< -► **tV 

— SUm 
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S * = ^ C «•£ _ ,<V* 0~ _ S’* •S/5- C ' ) 

_^ 0 t4^v- $(*■£ - ve W-* 


The ODE for the kth torsion node of the nth blade Is obtained by operating 
with Sr WK $*C***'>£**' s where 2% 1 8 the elactic torsion node 

shape. It is most convenient to apply this operator at this points 



&r — 

-^ A \*U.t-'otr < r-^^<.S ) <SUa J r' 



l.rf-n.tV'- s«. »*-**'>■*•* 1 

+ -VSV _ ( *.4 -».<•) i, J M 


and we shall use the notations 


y v - f (r-pa^ 


a) shaft motions 


y fF > *&* - (j'r.S 

+ 2 SI ( 

-f- ^ U 
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b) relative accelerate cnt 

r 1 ^ a*. 




where 


"♦* Ocb 

+ ( $’_Nt Vi f'~ a * s ) - (- + ?&"^b S ) 

+ Si <'.iu.iV ~SLf 
-*- L'** * + ^j.'x 1 © -ar 

— S' iL •(SiU-h^V r -* f >^©°~ ar 

J f PA> 

— Sr U© ^ ^ 

7 ’FA J 

s ^ •a *'*? 1 <&** — section pitch moment of Inertia, about SA . 


c) centrifugal acceleration: neglecting a number of terms due to blade 

torsion and pitch (of the same order as the 
propeller moment) , compared to the structural 
stiffening; there follows: 






- c 

+ E 9 (*+?•-**'<>'> &J- 

• ( x + 

t' > C *«%«r 7^ 

, C + "I 

'<>* r(w.-e-t^y*. Sr 
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elastic! usintf T = -St 1 . 

_ L -► 4 -n. 1 EX^©^ 3" 

d r 

+ Co^ tj si 1 £ >' 

+ - Ex tl £ } .(IotI-*„'4 ) ] 

aerodynanici the moment about the EA at r in 


thus 


a. . 5 ‘ «.*, >, . i 5 

S* j> . 

J?* = — — - 3*1 * ^ 

^L'* » **«•***, - 1 «- 

if 

« - ^ _ (*,-e - »«r£ v • ^ 

4 - *■ - - £> k * ■ + if A 


where 
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Rigid pitch 


The equation of motion comes from Mp^ •+ . where 

t* ?k~ V* - <~$9k + Cwe«*-*.*OWO m T* 

and M Is the nonent about the FA, at r = 


Inertia: 


M 


x o & a ** 


r* t r ( + <*.+»Vt +(*•♦»)*■ 1,^ a—«u- 

“ )Cp* ) V -(y.-t+^.TiMr,, J 


.So 


M ^r - ^ S L-(v?-^M rpA J 


+(tyt-W* 


&w£l r 


and we shall use the notation! 


a) shaft motion t 
I* 9 h m 


( ^' r ^. v^fiur^ ria-U 
^ c' *Sl. r«~ ar) 
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b) relative acceleration: 

A^Fa * - 


^ •• — 1 V 

) . ( — %Fa ©<* U -4* *fa e>c fc« ) 


_ c 5^X0*- -or) • (- 

c ' • ( » •*£ + ^ ^ ^ 

-’ rpA 

— C r © v^s^r 

}*>A V 

__ c‘ y. • < a, h» + (r " r f»'>®‘ '- tu * 

)Tfa 

+ 


c) centrifugal acceleration: 


Me*. = -JI 


^ S tP)k ^* ^ fcr ') • ** *«* 6 * 

+ ( S* *• -t. 

f» 

J<>A 

_ c' *\ • r ^‘*^3'“*^ 

)r ?K 

4 . 0 * ^.(x^ + ^ + ***Vfcr 

-*• u.*e ) c ^* ( *•■*- I oa ^ ** 

4. Curf 4 

— (‘ ~V*r 

*** I "<>* ^ffcl 

L |-*>fc( ijfft +&!%') J 
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aerodynamics : 


moment about the FA at r,,, is 

FA 

r* r 1 . n ^ <■' rrr-rn^^- ° c»- 

FA 




v rv.A *w^ 


where 


V.A„ - Yo 


elastic; 

The aerodynamic and intertlal moments about the FA are reacted 

by noments due to the deformation of the control system, due to commanded 

pitch angle, and due to feedback (mechanical or kinematic) from the blade 

bending or gimbal motion. The restoring moment about the feathering axis 

on the blade is — fl : it is given by the control system flexibility, i.e, 

oon 

the elastic deformation in the control system tines the control 

system stiffness K . Hence; 

con 

= ®cc * ^ce» ( £ — ®cm ^ ^ 

The q^ are the bending degrees of freedom, so Kp are the pitch/ flap and 
pitch/ lag coupling, mechanical or kinematic feedback due to the control 
system and blade root geometry. Similarly, Kp is the pitch/flap coupling 
for the gimbal motion. For the rigid flap motion of the blade, this 
coupling is given by the ftj angle, such that Kp * tan Sj . For a 


rigid control system, K 


, the rigid pitch equation of motion 


reduces to 


« *< 




So p 0 becomes just the control input, and pitch/bending coupling. 
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Now we write the control system stiffness K in terms of the 

oo n 

nonrotating natural frequency of the rigid pitch motion of the blade, 60 o , 
Then: 

Mp Ag - = CSc 

*4- £ V*f ; Y* 

+ ^Gr ^ 


Force 


The net force of the mth blade on the hub is 

F * *>-*■* 

where F is the force due to the blade, at the hub. 


s. \ 

inertia: = S. S "** fit/- 

* S«-i-^V®VN 

a) shaft motion: 

F = C s'o” 0 ^) — ° + 2SZ C ^„ rv “ tLr 

+ ( ^ o rw.ar') Ct^a-Tfatft’Ju. 

b) relative acceleration: 

p = (-7*. Hi -bit pc') -t- (>'.•<+ *.tV'~ar 

c) coriois acceleration: 

£ s 2Ji So S 
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d) centrifugal acceleration: 


— rt x S 0 S C'ta'tb r 


= - Si 1 


-£ # t So C ~* P * + ^wOr] 

+ 3* cC r “°^ -t « c So — ^ % 

+ S' *V + wSLr- 


aerodynamics : 


*> - 


Sur 


Moment 


1st 


\ -* 

sL f « 

* So L F »*» + F i^» 

•+ (Pr — + S»- S, -V X, '(*oi+-?o"^ T) J^*" 

The net moment of the mth blade on the hub, about the gimbal point, 

--J Cv*0 — » 

M — M* — A\ t 


inertia: ^*x “ So - <"** fi-X“ 

a) shaft motion: 

— c sl r ~®^) 

+ 252 C Si c*%*»£Lr) T* 'fo* 4 * 9 * 

-f* ( So «‘ p1 **®0 W# 
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b) relative accelerations 

/vs = 


c s'o fv ~ (*FA + >?A 0<j 'Zq') 

^ 9 C X ~~ 8L«- ^ c Vs “*• "t* (ids') 

( io~t — X o*£ } **"• ^-T* 

— S' e (( xo+^OT? + fcoll') r-v~SL r 


'f* 


-e°cs^i- $ r (r-<W}c war 


c) centrifugal acceleration 


*( 3 t*(St»p'>'> = c« ft il- ** = — rr l 

S C--e*T*-r> 


r 


.So 

fNA 


= n% 


r 


So C — ^F* * 4 ’ r S* — Cr-f>^Sj ") rwiSLr 

+ (Sl rlv ^ £Lr ) 

4. • (x#^ + *♦' ) r-v^SL^ 

4. (' © T& * (fc.t-x.tf -*x^} c»+&r 

- g° hMj ( slp A Cr-ci*)rm«i*) 
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d) coriolis acceleration: 

r x a s 2 Si. £ *ja X r ( r H'i -f -X^^ ^ 

+ ^»*vr "3 


cx/^x <i 


» - U (-x?a 4 4 -"ua 

4 - ( to't - y-oii - x*lt } 


,5a 


m - 2 SL 




-^fA ■*«■ it -(C-Cfy^Jj 


w&r 


4^6 |*-. 5 ^ r « 0 r 

j +^-(^-Xo^) (-r^fr-Cp^V^ 

4 ^ 0 T?l|'(*rf-««SV (*$*-*-<>* i)V^ 

4 S'o d V ~ fta* 

— (k G [-iPA ^ta.(»of-Xo^-X t l> 

+ fil ~ C^-^pX^ix J c** SLc 


aerodynamicr : 


M. - 


s 


' J c 
r x r, 

9 


ar 


a Vo ** “ R» ^ *) SLT 


-w- 


■• 1 
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Gimbal 


The equations of motion for the gimbal degrees of freedom are 
obtained from the ~C* t -$i components of M = £ M } thus: 


*‘HS + "i 


where is the spring and damper moment at the gimbal, reacting the 
rotor applies' •noment. The gimbal spring and damper are assumed to be 
in the nonrotating frame. Hence: 

= ( *<k “♦* 

— ^5 ( K<* (ic i»C. *4> 

_a 

Taking the T? 4 and components of M, the gimbal equations of motion 

are: 

-t* (bac (boc “ © 

+ lt(i (&G5 * O 


We shall write the gimbal hub spring and damper as: 

V.& a |1.^ Z ( 

C.& « ^ aco^x C(^ 


f * x 

where X* s > D r <*»&r 

flap notion. 


and ^ ir. tho natural frequency of the gimbal 
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Tlodal Equations 


Bending 

Consider the equilibrium of the elastic, inertial, and centrifugal 
bending moments, rrom the above analyrir, these terms give the 
homogeneous equation for bending of the blade: 

-*■ exx^it') y* 

_ jx x L sV J 

__ vm _st, • (to"i 

^.v>a C %o~d — ^ ~ ® 

This equation nay be solved by the method of separation of variables . 

Writing , % „ J\ — > . N \3t 

) =• v^(0 e. 

it be cones 

Cex^V”- sf CSr^®5 y ] ^ 

This is the modal equation for cour'n' 1 Hao/lag bending of the rotating 

o 

blade. It is an ordinary differential equation for the mo’e shape s 

this mode may be interpreted as the free vibration of the rotating beam 
at natural frequency ^ 

This nodal equation, with the appropriate boundary conditions 
for a cantilever or hinged blade, is a proper Sturm-Liouville eigenvalue 
problem. It follows that there exists a series of eigensolutlons 
of this equation, with corresponding eigenvalues . The eigensolutions 

— modes — are orthogonal with weighting function m t so if i ^ k , 

-I _i o 

5 0 t a-c *= © 

These modes form a complete series, so it is possible to expand the rotor 
blade bending as a series in the modest 

L - 1 

- 1 * 6 - 



1 


We shall normalise the bending modes to unit amplitude (nondlmensional) 
at the tip: 1 =1. 

Torsion 

Consider the homogeneous equation for the elastic torsion motion 
of the nonrotating blade; l.e. the balance of structural and inertial 
torsion moments, which from the above analysis ist 

— ©C, V* ■+ 31© ©e *0 

We could consider the equation for the torsion motion of the rotating 
blade, i.e. including centrifugal forces and some additional structural 
torsion moments. For the usual torsion stiffness of rotor blades these 
terms have little effect however, and the nonrotating torsion modes are 
an accurate representation of the blade motion. Solving this equation 
by separation of variables, write 6^, f SOS 

-CM j'V- X0 w l j *© 

This equation is a proper Sturm-Liouville eigenvalue problem, 
so it follows that there exists a series of eigensolutions , and 

corresponding eigenvalues (k * 1 ..,*•) • The modes are orthogonal 

with weighting function 1^ , so if i^ k 

The modes form a complete set, so the elastic torsion of the blade may 
be expanded as a series in the modest 

These modes are the free vibration shape of the nonrotating blade in 
torsion, at natural frequency • We shall normalize the torsion modes 

to unity at the tip, j^i) - 1 . 
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Expansion in nodes 

The bending and torsion motion of the blade is now expanded as 
series in the natural modes. By this means the partial differential 
equations for the motion (in r and t) are converted to ordinary differential 
equations (in t) for the degrees of freedom. 


For the blace bending we write 




where are the rotating, coupled bending modes defined above. These 

modes are orthogonal, and satisfy the modal equation given above. The 
q^ are the degrees of freedom for the bending motion of the blade. It is 
assumed (for the inertial terms) that the trim bending deflection is 
steady, independent of tlmei and when the substitution for the modal 
expansion is made, the subscript "trim" will be dropped, as that is all 
that will be meant by then. 


For the blade elastic torsion we write 

0C, — f ^ 1 1 ( r ") 

where are the nonrotating elastic torsion modes. These modes are 

orthogonal, and satisfy the modal equation given above. The p^ ( i ^ l) 
are the degrees of freedom for the elastic torsion motion of the blade. 


We also have the rigid pitch degree of freedom 

p 0 3 C6 # -e c ^ 4 

which is the total rigid pitch motion of the blade. Since it is rigid pitch, 
rotation about the FA, it has mode shape Si* Thus the total pitch 

perturbation of the blade Is as a series i 

IwO 
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For the blade pitch © then, the mean plus the perturbation is 

© — + © 2 C @t*U + d'tw') "♦* ^ 

0 

The subscript "m" on the trim pitch angle will be dropped when the 
substitution for the nodal expansion ir made, since that is all that 
will be meant by © then. 


ftonrotatlng Frame 


The equations of motion and the hub forces and moments are 
in the rotating frame yet. To get to the nonrotating frame, we introduce 
a coordinate transformation of the Fourier type i i.e., introduce the new 


degrees of freedom: 
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where (i© is the coning mode 5 \ the tip path plane tilt coordinates 1 

and is the reactionless flap mode — for the out of plane bending 

of the blade. Thent 

* £ C (^*«- -f (i*j -f. 

where the summation over n goes from 1 to (N-l)/2 for N oddj and from 1 to 
(N-2)/2 for N evenj the degree of freedom appears only if N is even. 

The quantities ,£*«., £* 4 , and p| are degrees of freedom, 

i.e. functions of time, just as the quantities are. These degrees of 

freedom describe the rotor motion as seen in the nonrotating frame, while 
the mP m ' > describe the notion in the rotating frame. 


] 


i 



This coordinate transform must he accompanied by a conversion of 

( rn ) 

the equations of motion for q from the rotating to the nonrotating frame. 
This is accomplished by operating on the equations of motion with the 
following summation oper •• rs : 


J. 

N 




N 




i 

) 7* 




Peference h gives more details of this transformation. 


Similarly, the degrees of freedom for the blade pitch ar.’ gimbal 
motion are transformed to the nonrotating frame. The corresponding degrees 
of freedom for the rotating and nonrotating frames are: 


rotating nonrotating 




i 

V s * >»<■>** 


Cl! 

p; 



£<Sr ^ ) $ aS ^ 

When the transformation of the equations and degrees of freedom 
is accomplished, there is a decoupling of the inertial and structural 
terms as follows (for N9- 3)s 

a) 0, 1C, IS degrees of freedom; . and ; 

and the rotor shaft notion. 

b) 2C, 2S, ... , ne, ns, N/2 degrees of freedom (as present). 
The first set couples with the fixed system motion. The latter set is just 
internal rotor motion. For N * 3» the first set is the conplete description 
of the notion of course. Nonaxial flow aerodynamics couples all the rotor 
degrees of freedom and shaft motion; i.e. the two sets above are coupled 

for helicopter forward flight or conversion mode operation. For axial flow — 
hover or pro pro tor airplane mode cruise operation — the aerodynamic terms 
decouple also. 

We shall assume here that the rotor has three or more blades, Nfc 3. 
For h • 2, there are periodic coefficients even in the inertia terms, so 
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that Is a special case. For the case of periodic coefficients in the 
aerodynamics, i.e. helicopter forward flight or conversion mode flight, 
it is necessary to specify N; we shall take N = R for that case. (The 
periodic coefficients depend on For the case of axial flow, or for 

the constant coefficient approximation for the nonaxial flow case, the 
equations obtained will be valid for all N greater than or equal to 3. 

Reference 4 discusses these points further. 


Equations of notion/ Hub Forces and Moments 

The elements are available now to obtain the equations of motion 
for the blade bending and torsion modes, in the rotating frame; and the 
forces and moments acting on the hub due to the mth blade. The steps 
required are: 

a) Substitute for the expansion of the bending and torsion 
motion as a series in the modes. 

b) Use the appropriate modal equation to introduce the 

mode natural frequency into the bending or torsion 
equation, replacing the structural stiffness terms 
(and for bending also some of the centrifugal stiffness 
terms). ^ ^ 

c) For the bending equation, operate with 

to obtain the ordinary differential equation for the 
kth mode of the mth blade (the q^ equation). 

d) For the torsion equation, operate with ^ 

to obtain the ordinary differential equation fo the 
kth mode of the mth blade (the equation). 

The result is the equations of motion and hub forces in the rotat!ng 
frame. The transformation to the nonrotating frame involves the following 
steps: 
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a) Operate on the hub force and moment with JL v ' . . . ^ ; 

i.e. sum over all N blader. to obtain the total force 
and moment on the hub. 

b) Find the components of the force 

and moment In the nonrotating frame (the S system). 

\ I> . * ■> 

c; '/Jrite the shaft motion a Q , U. , and uo in terms 

of the j^S components In the nonrotating 

frame (the S system). 

d) Apply the Fourier coordinate transform to the equations 
of motion and rotor degrees of freedom! operate on the 
equations for bending and torsion with 

-4 t . <■-•- 3 . k £ *— "> 006 ^ ^ ^ , ••• 

N ye * ^ 

to obtain the nonrotating equations of motion (0, 1C, 

IS, etc.). N ^ 3 is assumed for this transformation. 

The transformation of the equations to the nonrotating frame will be delayed 
however, so the rotating modal equations may be presented first. 


We arid at this point structural damping terms, modelled as equivalent 

viscous damping} the structural damping parameter is g (which may be different 

s 

for each degree of freedom), equal to twice the equivalent damping ratio. 


Names are given to all the Inertial constants now. The equations 
of motion, hub forces and moments, and inertia constants are also normalized 
at this point. The inertia constants are divided by the characteristic 
inertia 1^ ■ ^ r^m dr, and we introduce the blade Lock number ^ acH^/l 

This normalization of the inertia constants is denoted by a superscript *. 

The rotating equations of motion are divided by the hub forces and moments 
are divided by (N/2)^ for M x , M y , H, and Y, and by NI b for Q and T. The 
result is that the forces and moments are obtained in coefficient form. 

More details of this normalization procedure are given in reference 4. 


b # 
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Equations 

The resulting hub forces, hub moments, glmbal equations, and 
equations of motion for coupled flap/lag bending and for elastic torsion/ 
rigid pitch of the rotating blade are as follows. 
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Aerodynamic forces 
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and then for the nonrotating equations i 
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Inertia constants 
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Nonrotating frame equations 

The equations of notin’! for the rotor in the nonrotating frame, 
i,e. after application of the Fourier coordinate transformation, are 

+ A, Ao*fc -*• A t St -► A» o» •* Ao<* ** Svg + FA*. 


and the hub forces and moments 

• V •« V % A/ _ 

F * + C, ** -+■ Co%£ -I- C* oi +c, €>/-»• <^ex + r, 


where the rotor degrees of freedom ( x_, ) , shaft notion (a* ), rotor blade 
/ ^ \ 1_l -A 

pitch Input ( ), and the hub forces and moments ( ? ) arcs 
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The matrices of the coefficients , and the aerodynamic forcing vectors, follow. 
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' YY\ 'ICS 


In this section, tno urodyna r.\r forces arts monento op the r^tor 
blade are derived. We shall consider the general n:.e of high or low 
inflow, and axial or nonaxial t low . The aerodynamic terms in V'C rote »r 
equations of motion and the huh forces i nd moments are obtained for 
three cases: axial flow (hove^ nr high inflow cruise); nonaxin] flow with 

periodic coeixicients (helicopter forward flight, or conversion mode 
flight), and a constant coefficient approximation for nonaxial flow. 

The principle assumptions in the aerodynamic analysis are: 

reverse flow is neglected (good to an advance ratio of about 0 ,d or O.S 

which is sufficient for the tilting prourotor aircraft); the wing wake 

(near field an! far field) effect on the rotor, and other wing/rotor 

interferences are neglected; the unsteady rotor wake effects are neglected; 

the virtual mass aerodynamic forces and moments are neglected; the order 

c (rotor chord) terms in the aerodynamic lift expression are neglected ; 
a 

the order c^ terms in the aerodynamic moment expression are neglected; 
and only first order velocity terms are retained. The derivation and. 
notation are an extension of that in reference 4. 

Section Aerodynamic Forces 

A hub plane reference frame is used for the aerodynamic forces. 

All forces and velocities are resolved in the hub plane then, i.e. in 
the B system. The hub plane reference frame is fixed with respect to the 
shaft, hence it is tilted and displaced by the shaft motion. Figure 10 
illustrates the forces and velocities of the blade section aerodynamics. 

The velocity of the air seen by the blade, the pitch angle, and the angle 
of attack are: 

§ = blade pitch, measured from the reference plane 

u T , u^, Up ~ air velocity seen by the blade, resolved with 

1 respect to the reference plane; u^ is in the hub 

plane, positive in the blade drag directions 
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where 


is in the hub plane, positive radially 
outward along the blade; and Up is normal 
to the hub plane, positive down through 
the rotor disk. 


U = resultant air velocity in the plane of the section. 
= induced angle 
OC = section angle of attack 

, X T * 

^ ss Tsx ^ ' lAf/^-r 

OL =z © — & 


The aerodynamic forces and moment on the section, at the EA , are: 


L, D = aerodynamic lift and drag forces on the section, normal 
and parallel to the resultant velocity V 

V ^ , F = section L and D (total aerodynamic force on the section) 
resolved with respect to the hub plane, normal to and in 
the plane of the rotor 


r radial drag force on the blade, in the plane of the disk, 
positive outward (the same direction as nositive u^); 
the radial forces due to the tilt of i and F have been 
considered separately. 

- section aerodynamic moment about the EA, positive nose up. 


Aerodynamic forces — wind axes 

The section lift and drag are 

L. s' 

*5 Wee* 

where U ■ resultant velocity at the section 

> ^ = air density 
c “ chord of blade 

The air density is dropped at this point, in the process of making the 
quantities dimensionless with o , , and H. The section lift 
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and drag coefficients, , , ) wid = c^( *<r , il) are funct'or: 

of the section angle of attack and Mach number: 

c* s a - 4> « g> -Tcur’ 

M — h*-r)P U 

where M ^jp i- the tip Mach number, the rotor tip speed jfl : ' divi *ed by 
tne speed of sound. The dependence of and c^. on other quantities, such 
as the local vaw angle or unsteady angle of attack changes, is neglected. 

The radial force, due to the radial drag, is 

F r = ^ 

r UA ^ 

The radial drag force is derived assuming that the viscous drag force on the 
section has the same sweep angle as the local section velocity. The 
moment about the 3A is 

^ Sr — >a L - *♦” m AC ♦ M 

m — >1* C-C.fi *+* lA l C. 1 ^ 

where x, = distance aerodynamic center (AC) behind EA 

A 

c - section moment about the AC, positive nose up. 
m 

ac 

- unsteady aerodynamic moment. 

lor the section aerodynamic moment it is necessary to include the unsteady 
aerodynamic terms, which from thin airfoil theory are 

. £*r cvavi+t 

® c ji u 

•b ( li + uai* ^ C \ b ■ t ^ 

where w = mean upwash along the blade chord, i.e, normal blade section 

= u^slnft - Upcos© 

B = , basically the pitch rate ^ 

V - i^cosd + UpSin® 


If 
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Hence In the aerodynamic model we have neglected the following 
effects: reverse flew; shed wake aerodynamic interference lift 

deficiency function set to unity); term:' in L order c .and above; term: 
in M order c and abovv ; virtual mass terms in the unsteady aerodynamic 
moment. 


Aerodynamic forces -- hub plane axes 

With respect to the hub plane then 

F-2 - LoO-S^ ~ ts 4 


Lu t - 
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U 

Substituting for L and D, and dividing bv ac (where a is the two-d imensional 
section lift curve slope, and e the section chord; which enter the Lock 
number Y also), we obtain: 
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The net rotor aerodynamic forces are obtained by integration of the section 
forces over the span of the blade, an then summation over all N blades. 


Perturbation forces 

Bach component of the velocity seen by the blade has a trim term, 
due to operation of the rotor in its trim equilibrium state; and a perturbation 
term due to the perturbed motion of the system. The latter is due to the 
system degrees of freedom, and is assumed to be small in obtaining the linear 
differential equations describing the dynamics. We shall write the blade 
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pitch and section velocities as trim plus perturbation terms: 

& *o e -4- £e 

l«T U*r -f £v*t 

Uf aa(> Up -+■ «W 

lAft =l> •* 


then there follows the perturbations of 04 ,U, and M: 

So = 6 © — UT &*p - U^U T 
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u x 

iat 6 i*t 

C* 
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and of the aerodynamic coeffients 

S<* -V So* -4- CJ M 

(and similarly for c and c.). The perturbations of the section aerodynamic 
forces may then be obtained by carrying out the differential operation on 
the expressions above for F , F , F , and M , using the above results 
to express the perturbations in terms of £& iSu^, Sup, and Su^. 

The coefficients of the perturbation quantities are then evaluated at the 
trim state. The results are » 

*sl~ CUU-ryg - LAI** 
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Velocity of the Blade 

Now we obtain the velocity of the air seen by the blade section. 

There is the trim velocity, composed of the forward speed, rotor rotation, 
and rotor induced velocity; and the perturbation velocities, due to the 
rotor degrees of freedom and the shaft motion, and due to the aerodynamic 
gust velocity. 

The rotor is rotating at constant speed -5~1- . The steady velocity 
of the rotor with respect to he air, is described by (figure 11 )i 

V - trim velocity of the rotor in inertial axes, 
in the rotor x-z plane. 

angle of attack (undisturbed) of the rotor hub plane 
with respect to V, positive for disk tilt forward (for 
V down through the disk); this is the shaft angle. 


There are then the following cases: * 90° for cruise (high Inflow 

axial flight)} small for helicopter forward flight; c* w p large 

but less than 90° for conversion mode; ?nd V « 0 is the hover case. 

The rotor induced velocity Is v, due to the thrust T (figure 11); v is 
assumed to be normal to the hub plane, and uniform over the disk. Now the 
rotor advance ratio f* and Inflow ratio \ are defined} 


A “ 


N/ ccs tw *« 


V fvw ©'vP + v 

The cases are then: for hover /+ m 0 and \ small: for helicopter forward 

flight 0 and V small; for conversion mode flight 0 and 

\ order 1; and for cruise flight f* - 0 and V, order 1. 


For the rotor Induced velocity wc use the Glauert result; 


\ b +* 


For high speed ( ->■> ^C^StR) 4 * or about V/TL R > 0.15) in inflow ratio 

Is approximately 


jtk + z v/i.fc 


The induced velocity is thus quite small, v/V « 1, for typical pro pro tor 
cruise and conversion mode operation. The induced velocity is not generally 
an Important factor in p ro pr o tor aerodynamics at high inflow; hence the 
assumption of uniform induced inflow is acceptable for an investigation of the 
pro pro tor aeroolaatic behavior. (See reference 4.) The mutual aerodynamic 
interference of the rotors is neglected. 


The trim velocity V le steady, at an angle ©<hP to the rotor hub 
plane. The uniform induced velocity v is normal to the hub plane. The 
advance ratio and Inflow ratio, and X , are the nondinenelonal 
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components parallel and normal to +!,e hub plane. In body axes, V would 
be fixed in the reference frame, and would tilt with it. Here an inertial 
frame (the S system ) is used however, :•< it follows that tilt of the rotor 
by the shaft notion gives a s -11 change in the direction of V as seen in 
the reference frame. 


The shaft motion consists of snail linear and angular velocity, 
with components defined in the nonrotatirg frame: 


an* = + y* Tl* 

*We — ^ + •'a 



The aerodynamic gust velocity has components u^ t v^., and w^, (longitudinal, 
lateral, and vertical) defined with respect to the body or earth axes 
(figure 11) j the6e components are the velocity seen by the aircraft, and 
are assumed to be small compared to -fl-H . The gust components are 
defined with respect to V, i.e. from the disk plane, so that with 

V usually horizontal (level flight) w Q and u^, arc always the vertical 
and longitudinal components with respect to the flight path. The gust 
components are normalized by dividing by J1.R, not by V as is often the 
convention for airplane analyses. The aerodynamic gust is assumed to be 
uniform throughout space. 

Trim terms 

The result for the trim velocity terms ist 

Ur « «■ +/* ^ — /-uiV ^ 

L V* T* 

V + ty-ti-T + 
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Ufc - o©sV + ( VRA *♦* "j *+* C *V* * ( v { ; ^ ^ 

—v ( &*< t -\vA x + "** (* 0 -^ 

+a>^ ( - s-vT**^ } 

CK/m^ 

0 — ©c®\' -V- 9 tv» 1- ©C*y. ~ ^<r — LY-fl 

Wiere ©o^t. is the input cyclic pitch required to trim the rotor. 

“or the trin velocity, the blade bend! and girabal notion 5r pericdic. 

’.•or axial fl : rht, /V* = 0, the trim velocities are constant; for nonaxlal 
flow, y* ^ 0, these velocities are periodic in , due to the rotation 

of the blade with respect to the rotor forward velocity. 

Perturbation terns 

The result for the perturbations of the velocity components and 
the blade pitch, due to the rotor and shaft motion and the aerodynamic 
gust, is then: 


&\*T — CVo '* •+ -*■ C-041 V~. 

•4* ( — Xv* + 4 S*** V*** 

4* txti Hw ( ©f g, 4 y$ ^ 

■4 r -4 ^4 ^ 

-4 L (£4 • Y ") 

4- y* £ «v* 



£{a* =* 


T" yMt«S *V*« C 

■+ C ( “h O/y i^HU* ” C^tj e»4 

+ £. <v 

'+"yucoi*V— £ 'V (■ * I 1 ^ 

<§Ufc = — ( V. + \^ u -4- ^ Vv-» 

■4* ( Vo<*j — >v» -+” U<> t*A<V H ^ -fV9(; coaV*** 

— 4 ^ 

■+■ £ £ lcj| • (~2^\ — f — yM f>M H*» *|'. ^ ^ 3 

= a = £7 pi *i 

C*r MwS 

vB = (ut cj +*&-*- ia^ fw e> ( d + (J<* *+ £^* / 

u>4u fc io r * Z pt J* (uTM*e4.K^e) 

■4*^(1 £plC JiT < u T ***0-»'WfM?)4’ 2|: «^0') 

— <l<* (2u*t»*er> 

•+- £o (/*s**V«* c^d©^ 

-+ <©<j -v C^UfcSM.®^ 

— C®*4 -fSfj^ C *** 

— 2u* £<i; 

-h £ Y* 0 % “*~ ^*T ~ u * *‘T ' 
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Rotor Aerodynamic Forces — jtojtetln£^_Blade 

With now the expansions for the section forces and monent in terns 
of thd^elocity perturbations, ami the velocity in terns of the notion 
of the rotor and shaft, we nay obtain the perturbations of the aerodynamic 
forces on the blade. These are the blade forces expanded as linear combination^ 
of the degrees of freedom. Giving names to the aerodynamic coefficients 
at this point, the vesults for the required aerodynamic forces on the rotating 
blade are as follows. 

Bending: 

"*’ V ° V ^“ 4 

4* O 

+ NV^ k \ ( 4k —/* Ot^ -4 U* “ w * ***° / »^ 

4> K*»*Vi»> “• ®<vj 

■*- £ V 

4 - £ 

*+ £ 
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radial force: 


So si - si 

~ ^ — (Vo/jj 4(^u 

+ (V<*<*3 —Xu -4- ^c^e4 H *4-*0*fc~c* H ^«»* , K^'] 
+ 2r C (*<** + *£* + ''* ^ c -^ 3 ^' 

-4 c >OAj — Xu -4* + A** SWV*."J 

-4- (£/* + H0 -4- &£ (cx* 4-^ 

4- fcy ( JJu —A “♦• **• t>**o>*t — Wfr c *4°^iP >S ) 

*+ ( (i* 4- cv* K*HU — <^vj + fcp ^c, 

4- £ <^5 V -+• r* 

4- t&f. {* 

Torsion/pitch: 

sL w**^ - c c &* + aj*> -v ^ 


— 4- vy, 4V^CW 

4- (V*^ -X*» ♦ M&co*cu w ^ 4 wifi, 

+ ^i') -4* C +H^ 

4* M^vlV L4v — /* <V%J -4* U*fuo(^ — 

•fo»* 

+ l J*e*v ■+ ^ M f*v T 

"♦" C f i ■* £ P : 
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Hub forces and moments: similar to the bending case, but with notation 

integrand notation 


flap moment 
torque 

blade drag forcr 
thrust 


rF 

2 

rF 


h 

Q 

H 

T 


Aerodynamic coefficients 

Applying the results for the expansion of the aerodynamic forces, 
and the expansion of the velocities, the aerodynamic coefficients may be 
evaluated . These coefficients of the degrees of freedom in the aerodynam* c 
forces are constant for axial flow, the J* = 0 case* For the general nonaxial 
flow case, /t* ^ 0, the coefficients are however periodic functions of . 
The results follow. 


Bending: 



So* 

- 

si 

*%i - 

Vo 


r*- 

* 

si 

- 

si 


A * 4 

\ 

r**\*v * 

s. 


\ cc Cp- ^ 

tft ~ &T 






/acAa'VU Q So SlS 

Slt 


\ 



Flap moment: 


Ny 

- So p ir <-«U- 


'•i* 

» s: f*- 


M S 



A*v 

• Si 



- Si Fir <■'** 


= r u.f„ (A X 


M r- 

* sic«w4- 

‘T* <-* ♦ Y ) <“SLr 

N^i 

So 

^t"^*^** x?% * Y* v } rfi>r 


= si F 4e L 

cSU" 


Other hub forces and moments: similar to flap moment, 

coefficient 

flap moment I* 

torque 

blade drag force H 

thrust T 

Radial force: 

S* = S'. ** ^ 

^ = si -Hr U.~ »t n &r 

<4 = C CSU 

(ty » S*» t*»y — F*p (Si —4* +&*(*»T!+ , 4 o^ ^ 1 

» Si tP«r +ii-(u-t+ut'>’ , n <■&*■ 


with 

Integrand 

rF 

z 

rF x 

F 


-as- 



) 


= — So ^ &- i " 

- C + .(we+Wi^^l^-Y <Slr 

■* So £^«> -F iP Ur£L+T^*(*<»^+%©"£V 4] x**^ &r 

•+S. 

-v So^<“r L^* (*f i — r »Ji — /•t^Hu y v4*y 3^" 

“ Si ll T* -^iSU 

Kp*. — ^ 0 -" ( 5 i“^j>'*''^ * (voT?-f ^^^;£U“ > 

Torsion/pitch: 

r 1 r , 

Mfn/* s jfp^ [ |»^r — (^xt "3 8lT 

= $lp A t Kk-p - ( SSc T fa + P * T 

Srp A 1 x J 2 (l+ 4 3 2^1 v~<? 

M f*£ = ^.c^H’U 

■+• 4 r* ^ f£ ( ,+ j »rV^^ s*.©^ 
**»> * ${** t Jh Kf — -i- P^Tc* V^“j£Lr 

- $c*a It*. *^e “ (^»f*6i T rStr" 

+ ic^ixu (*♦ 4^^ 2 ucc*a© ©r 

*hfc* 

— S\** Uji (uT^»+krv*-eSOi-t 
— Sfy A 3 a fa 0* 4 ^3 y*fc~Hu c%*9 Sir 
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+ s)j. A +^if^»')-\f.~l 1 **’T ft/ 

+ C* UJiO+'l 1** *•$’*'' 


t*f*V = 




^-V- \ W C-U^T -<F,rt»+fST^>JV"i^'f ^ 

+ S^» U***f - < F »^* +f *^V lf *T ^ 

- SW '“W *** a , 

-S'cs. u& *««-•♦ >«♦*■• ®Vi+t % H'T ^ 

- s' r(>A \v. £ O+N^Jf*- 1 * *T " (r^' ) ’" e 'f^ Su ' 


s 5 *^ tj„>«e- (F« e f»+ F *et»VW 31 i ^ 


Htfx - -Sr^UV € < “— 

a - (Uf +Utt|+Ajt|^( r "^^ 

+«.<-*.'*') I^+ 

Wl» = 


- 90 - 




1 


Rotor Aerodynamic Forces — Nonrotating Frame 

The aerodynamic forcing functions for the rotor equations of 
motion in the nonrotating frame, and the hub forces and moments are 
now required . These are obtained by summing the blade rotating forces 
(given above) over all N blades. The Courier coordinate transform of 
the rotor degrees of freedom is introduced as required. 


Axial Flow 

First consider the case of axial flow, /* = 0; for either 
high Inflow ratio \ (order 1, i.e. proprotor cruise flight), or low 
inflow (small \ , i.e. hover in helicopter mode). In this case the 

aerodynamic coefficients in the blade forces are constant, independent 
of hl . The coefficients are also independent of m (the blade index) 
then, so the summation over N blades operates only on the blade degrees 
of freedom and shaft motion variables. The result for the required 
aerodynamic forces, in matrix form, is 


f v - 

* A, Kfc 4* A,*fc + A» o( + As 

F<xju» » c i + C*** ♦ bfc \ 


where the rotor degrees of freedom ( Xp ), shaft motion ( o< )» and 
aerodynamic gust input ( g ) are* 


1 



9 h . 

(&«« 


Of s 


~Xv» 


— •* 

3* 

-4 

Uo- 

2w 

cv* 

L^J 

y 

Vo 

Wo 


These coefficients simply add to the inertial coefficients already derived, 
to complete the equations of motion. The matrices of the aerodynamic 
coefficients follow. 
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Nonaxlal flow 

Consider now the C3.se of nonaxial flow, /* > 0, This case 

includes helicopter mode forward flight, and conversion mode flight for 
the lilting proprotor aircraft. The aerodynamic coefficients are then 
periodic functions of . Hence the equations of motion for the 

system have periodic coefficients, due to the periodically varying 
aerodynamics of the edgewise moving rotor. 

One can express the aerodynamic coefficients as Fourier series, 
and then obtain the coefficients of the nonrotating equations of motion 
in terms of these harmonics. For the general rotor considered here, 
it would be necessary to evalaate the harmonics of the aerodynamic coefficients 
numerically, however. It is simplest therefore to just sum (numerically) 
the coefficients over m * 1...N as is required in finding the nonrotating 
equations of motion and the net hub forces and moments. The nonrotating 
coordinates for the rotor motion (Fourier coordinate transformation) are 
also introduced. 

For the periodic coefficient case, it is necessary to specify the 
number of blades N, since the periodic coefficients depend on Ns also, 
the periodic coefficients couple all the rotor nonrotating degrees of freedom, 
so more than the 0, 1C, and IS variables are involved with the shaft 
motion (if N > 3)» We shall consider only the case N - 3» then the 
0, 1C, and IS degrees of freedom are the complete set, even for the 
periodic coefficient case. The period of the equations in the nonrotating 
frame is AH' * 2w/N. 

Again we write the aerodynamic forces in matrix form, as 

iw ^ 

— m + k\ 

fWe ■» c » i* + + & •* + 2s a* 4* 
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where now the coefficients A, D, C, and D are periodic functions of 
(period 2 tt/N). The matrices of the aerodynamic coefficients follow. 
The notation 

C = cos 
S - sin MU. 

is used (Mu s . /wV ). Note that each matrix is a summation over 
all N blades (N ■ 3 in this case). 
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Constant Coefficient Approximation 

Finally, we consider a constant coefficient approximation for the 
nonaxial flow case. This approximation uses the mean values of the 
periodic coefficients of the differential equations. A constant coefficient 
approximation is desirable (if it is demonstrated to be accurate enough) 
because the calculation required for the analysis is considerably reduced 
compared to the periodic coefficient equations, and because the powerful 
techniques for analyzing time-invariant (constant coefficient) linear 
differential equations are applicable. It is only an approximation to the 
correct dynamics however; the accuracy of the approximation must be determined 
by comparison with the correct periodic coefficient solutions. 


To find the mean value of the coefficieats, we apply the operator 

$lr ' ) ©* v 

to the periodic coefficients given above. The result is terms of the 
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where M nc,ns are the harmonics of a Fourier series representation of the 
rotating blade aerodynamic coefficient fit 

M * r*<° -4- £ m < ***%»H' % * •• 

v»* » 


In the present case, these harmonics must be evaluated numerically. We 
evaluate M at J points, equally spaced around the azimuth t 


^ EM; 




where 

*»i « i** j i='" ' ar 

A M> x ^ 


The harmonics up to the second (n = 2) are required here. This Fourier 
interpolation formula requires then for good accuracy about J> 12 (a 30° 
azimuth increment). Using these expressions, the required harmonics aret 
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- £ M U 


4m 


it 


£ 


/. 
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C04V 
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2-Sw** j 
2. V**H*«ja*H > / 


**00 
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It follows then that the constant coefficient approximation is 
obtained from the periodic coefficient expressions by the simple 
transformation : 


K? 

W wa» 




The summation over N blades (n = 1...N, AV = 2W/N) for the periodic 
coefficient case is replaced by a summation over the rotor azimuth 
(j = 1...J, 2-lt/j) for the constant coefficient approximation. 

This is quite convenient, since the same procedure may be used to 
evaluate the coefficients for the two cases, with simply a change in 
the azimuth increment. The periodic coefficients must be evaluated 
throughout the period of *P - 0 to 2ir/N of course; while the constant 
coefficient approximation ( the mean values only) is evaluated only once. 

With the substitution j 3T j , the results 

given above for the periodic coefficient matrices are directly applicable 
to the constant coefficient approximation as well. 
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ROTOR TRIM 


The are two requirements in the dynamics analysis for the trim, 
equilibrium solution for the rotor blade motion and rotor performance i 
first, the trim bending deflection is required for the 

coefficients, particularly when the bla-’o torsion dynamics are involved} 
secondly, the evaluation of the aerodynamic coefficients requires the lift 
and drag loading of the rotor blade. The trim bending deflection is assumed 
to be independent of in the analysis, so the mean value must be used 

when r ' 9 0 ; for the aerodynamic coefficients, the periodic variation of 
the trim blade aerodynamics when 0 will be included however. 

The dynamics analysis (the evaluation of the coefficients of the equations 
of motion) must be preceded therefore by a preliminary calculation of the 
rotor equilibrium motion. The trim solution for the blade motion is 
periodic in the rotating frame for the general case of nonaxial flow i 
for - 0, axial flow, the blade notion is steady in the rotating frame. 

For the trim blade motion solution we shall consider only the bending 
and gimbal degrees of freedom. It is assumed that there is no shaft motion, 
gusts, rotational speed perturbation, or torsion/pitch motion (except 
cyclic control and any bending/ tors ion coupling) in the trim solution. 


The trim solution involves the numerical Integration of the 
differential equations of motion for a single blade in the rotating frame, 
until the blade motion converges to the desired periodic solution. The 
equations for the blade motion are obtained from the above analysis, and are 
for the bending and gimbal degrees of freedom 

* + * si 



£ c 
ar j 



-ii?- 



where the inertia constants are defined above, and the aerodynamic forces 
are evaluated using the trim velocity components (for which expressions 
aro given above). 


the integration of the blade notion converges to a periodic 
solution, the rotor performance may be evaluated, i.e. the mean aerodynamic 
forces and moments the rotor produces at the hub, particularly the rotor 
thrust and torque coefficients. The Fourier harmonics of the blade bending 
motion are also evaluated. From the zeroth harmonics oi the bending motion, 
the mean bending deflection of the blade may be evaluated. 

• or axial flow, Ja • 0, integration of the blade motion is not 
required; for the gimbal motion is zero (assuming no cyclic pitch input) 
and the equation for the blade bending modal deflection reduces to 

■=-v» V + V 



ELATE BENDING AND TORSION MOTES 


Cowled Lending modes of a rotating blade 

Equilibrium of the elastic, inertial, and centrifugal bending moments 
on the blade gives the differential equation for the coupled flap/lag 
bending of the rotating blade. Kor free vibration — the homogeneous 
equation (no forcing) with harmonic motion at the natural frequency ^ — 
we obtain the modal equation for bending of the bladet 

csz Y" — si 1 ( 5^ ^ — ***^*^*° 


B* - 
SL st 

Q 


where « bending deflection (mode shape) 

ex - m bending stiffness dyadic 

st s Ht| » rotor rotational speed 

^ m natural frequency of mode 

This is an eigenvalue problem, a differential equation in r for the 
mode shapes ^ and the natural frequencies ^ . The equation with 

the appropriate boundary conditions constitutes a proper Sturm-Liouville 
eigenvalue problem. It follows that the solution exists t a series of modes 
and corresponding natural frequencies ?; | where the modes are 
orthogonal with weight m, i.e. if i^tk then 


^ so 


and the frequencies satisfy the relation (an energy balance) t 


The modal equation will be solved by a Galexkin method. The mode 
shape is expanded as a finite series in the functions f^(r)t 

fd?i 
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We require that each of the satisfy the boundary conditions on ; 
then the sum automatically does. Since a finite series is required for 
computation, this is an approximate calculation; the functions f . should 
then be chosen so that at least the lower frequency modes can be well 
represented, for best numerical accuracy. Substituting this series in 
the differential equation and operating with 

S), V<--OAf 

reduces the problem (after integration by parts and an application of the 
boundary con ' tions) to a set of algegrMc equations for c - 

C K - -O' 1, & } c so 

where the coefficient matrices are 


B k *. — 

't 


Eigenvalues of the matrix 3"*A are the natural frequencies O 1 of 
the coupled bending vibration of the blade; and the corresponding eigenvectors 
c give the mode shape vg • As a final step, the modes are normalised 
to unity at the tipi - 1, 

A convenient set of functions for are the polynomials l') t 

“T* * to 3 TT 

(for a hinged blade f ^ ■ r is used). These polynomials satisfy the 
required boundary conditions, but are not orthogonal functions. 
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Torsion nodes of a nonrotatlrr blade 

Equilibrium of the elastic and Inertial torsion moments gives 
the modal equation 

can yy * j«o 

The modes are orthogonal with weight ; i.e. If 1 fc then 

C - o 

and the frequencies satisfy the relation 



These are nonrotating modes, no the solution is independent 
of SI or & . The equation is solver! by a Galerkin method* Writing 

CC! 


where the functions satisfy the boundary conditions on J , and opera ing 
with n<‘"^ke on the differential equation, produces a set of 

algebraic equations for “c - T « 

(A - 4*0^ * O 


where 


A *‘ “ $*»*. 


&*•. - 5c* = ^ 
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The eigenvalues of the matrix D _1 A give the natural frequencies of the 
torsion vibration, and the corresponding eigenvectors for c give the 
nodes. Finally, the torsion modes are normalized to unity at the tip, 

- 1 . 

A convenient set of functions to use for f^ Is the solution for 
the torsion nodes of a uniform beam: 

r 1 \ «'“*’•* -i 

-f. * a- C f— J 

These functions satisfy the boundary conditions, and will usually be 
close to the actual mode shapes. 
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SUPPORT EQUATIONS CF MOTION: CANTILEVER WIuG 


For the rotor support we consider a cantilever wing, with the rotor 
on a mast or pylon attached to the wing tip. Reference 4 discusses the 
cantilever wing as a representation of the tilting proprotor aircraft 
dynamics, and develops the equations of motion describing this support. 

The equations of motion for the wing, and the rotor motion produced by 
the wing are developed In reference 4; these results are adopted here 
with only two extensions: to arbitrary angle of attack of the rotor 

shaft with respect to the forward velocity*, and the Inclusion of a wing 
t railing-edge flap among the controls. 

Cantilever wing 

The cantilever wing and pylon geometry is shown in flg'.r'e 12. 

Ve consider a high aspect ratio, flexible wing, with the rotor on the 
tip. The wing is attached to an Immovable support with cantilever root 
restraint. A pylon with large mass and moment of inertia is rigidly attached 
to the wing tip. The rotor is mounted on the pylon with the hub forward 
of the wing EA, with mast height h. A general pylon angle Spis considered, 
from vertical in helicopter mode to horizontal in airplane mode. The wing 
motion consists of elastic bending, vertical and chordwlae, and elastic 
torsion. There is no motion of the pylon relative to the wing tip, so the 
wing tip motion is transmitted directly to the hub, and hub forces and 
momenta transmitted directly to the wing tip, through the mast of height h. 
The rotor and wing operate in a steady free stream of velocity V. The 
pylon (or asst, or rotor shaft) angle of attack £f» may be large, so it 
covers the entire range of tilting proprotor operation. The cases includes 
6 ? near 90° for helicopter mode: if between 0 and 90° for conversion 

node i * 0 for cruise node; and V * 0 is the case of hover flight. 

The wing angle of attack is , defined positive nose up: 

It is assiaied that 4e 4 la a small angle. The angle between the wing and 
the rotor shaft la then i It la this angle which determines 
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the transmission of motion and forces between the rotor and the wing. 

Recall that is the angle of the rotor di 6 k to the forward speed V 5 

here we use <5f> for the shaft angle of attack, hence <5, = 90°- . 

We also consider small sweep angle j (positive aft) and small dihedral 
angle iw, (positive up) of the wing. A major effect of Swj and is 

on the position of the effective elastic axis of the wing, hence on the 
effective mast height for the transmission of motion and forces between 
the rotor and the wing. The angles Jw, , <Su t , and 5-^ are removed from 
the orientation of the pylon and shaft at the wing tip. So the rotor shaft 
is in a vertical plane with no sweep or dihedral, parallel to V when 

= Oj ar.u then is the angle of attack of the shaft with respect 

to V, not with respect to the wing. 

The wing is assumed to have a straight spar line, which is the locus 
of the local iSA. The wing root is supported with cantilever restraint, and 
the rotor shaft is attached rigidly to the wing tip. The wing has no twist, 
constant chord c w , length y T from root to tip (semispan), with the distance 
y w measured from the root, along the wing spar. The shaft length (mast height) 
is h, the distance the rotor hub is forward of the wing tip EA. The wing 
spar is roughly perpendicular to V, with small wing sweep, dihedral, and 
angle of attack considered. The wing root is attached to a plane defined 
by the forward velocity V and the vertical; then the three rotation angles 
Sw, . . and define the orientation of the spar with respect to 

the free stream velocity. Next the pylon is rotated by iw, . , 

and with respect to the wing tip, to keep the shaft parallel to V ; 

finally the pylon is rotated by with respect to V , defining the 

orientation of the rotor. 

Swept wing6 are usually built with a center box structure in the 
fuslage, where the spars are unswept, and only the wing structure outside 
the fuslage has swept spars. The wing is restrained at several points where 
the wing box is tied to the fuselage structure, or in this case to the 
cantilever wing fixed support. There exists an affective elastic axis 
for vertical bending of the wing tipi some point on the shaft or its 
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extension where the application of a vertical force results in purely vertical 
displacement of the shaft, with no rotation in pitch. Without sweep 
this point would be just at the wing tip EA? but with sweep a force there 
will produce a pitch motion of the shaft also, hence the effective EA is 
some distance from the wing tip EA. The effective elastic axis for the 
tip lies between th^LCtual wing tip EA and the extension of the upswept 
spar line, the actual position depending on the degree of root restraint 
and sweep, arid other s+ructural details. Figure 13 illustrates the 
geometry involved. Reference 4 develops an elementary model for the 
wing bending and torsion including the shift of the effective EA due to 
sweep (and a similar effect due to dihedral), which is adopted here. 

The effective EA position is described by (figure 13)* 

h = mast height, distance hub forward wing tip EA. 

h_. = effective mast height, di tance hub forward 
effective EA. 

= distance hub below effective EA due to dihedral. 

Further discussion of this effect, including the estimation of the 
parameters involved, is given in reference 4. 


The aircraft has two contrarotating rotors, one on each wing tip. 
The direction of rotation of the rotor on the right wing (as in figure 12) 
may be either clockwise or counterclockwise. The influence of the rotor 
rotational direction is a few sign* in the equations of motion, reflecting 
how the rotor hub forces and moments excite the wing motion, and how 
the wing produces motion of the rotor shaft. As in reference 4, the 
notation -TL is used to carry this influence of the rotor rotation 
direction, where S~L takes only the values ds 1 1 





+ 1, rotor rotation clockwise on right 
wing, counterclockwise on left. 

- 1, rotor rotation counterclockwise e i 
right wing, clockwise on left. 
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Wing Motion 

The wing motion in ’escribed by elastic bending and torsion of the 
spars the pylon, and with it the rotor shaft, is rigidly attached to the 
wing tip. Elastic bending results In deflection of the wing spar with 
components both perpendicular to the wing surface (vertical or beam bending), 
and parallel to the wing surface (chordwlse bending). Vertical and chordwlse 
bending are 'efined with respect to theji tree Lion of the local principle 
axes of the section. There is no wing twist, so these principle axes are 
the same all along the span, but they are not vertical and horizontal nxes 
because of the wing sweet), dihedral, and angle of attack. We define (figure 
IT) the wing bending and torsion deflection as follows: 

z (y ) elastic bending vertical displacement of the 
w r ;xi r , normal to the wing surface, positive up. 

x (y ) = elastic bending chordwlse displacement of the 
M w span, in the plane of the wing, positive 
rearward . 

& (y ) ~ pitch change of local wing section, due to 
w w elastic to rs ion about the local EA, positive 

nose up. 

A modal description of the wing elastic deformation is used, and only 
the lowest frequency modes retained. Wc consider just three degrees of 
freedom for the ring: first mode vertical bending, chordwlse bending, and 

torsion. The degrees of freedom representing the wing notion are: 

q * wing vertical or be&mwise bending, positive 

W 1 upward? q = t /y_ at the wing tip. 

W 1 w 

a “ wing chordwise bending, positive rearward } 

* % m vV T at the win « ti P* 

? A w 

p ■ wing elastic torsion, positive nose upj 

W p “ ^ at the wing tip. 
w w 

Associated with these degrees of freedom are mode shapes, fo 1 " 

torsion, and for bending. These modes are normalized to 1 and 

to y T respectively at the wing tip (y y - y T ). 
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Front the results of reference 4, generalized to arbitrary pylon 
angle of attack , Sp , the rotor hub motion due to the wing degrees of 
freedom is: 
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where we have written 


G for cos( 

S for sin(£f — J>o^) 
•J for <•}**•* 

y for y T 

w 

for A Jo I 
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Equations of Motion 


r'rom reference 4, the equation of motion for the o , q , a nr 

W 1 M 2 

m degrees of freedom of the cantilever wine;, excited by the forces and 
moments at the rotor hub and by the wing aerodynamic forces, are: 
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The wing equations are normalized by dividing by (N/2)l^ # so th^rotor exciting 
forces are in helicopter coefficient form. The inertias are* 
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where m w is the wing mass per unit length; I^ w is the wing section moment 

of inertia in pitch; »p is the pylon mass (without the rotor); Ip and Ip 

x y 

are the pylon yaw and pitch moments of Inertia, without the rotor', about 

the wing tip effective EAf and z v is the distance the pylon CG (without 
. *ea 

the rotor) is ahead of the wing t3p effective EA. For the pro pro tor 

configuration, the pylon mass is so large that it dominates the wing 

inertias. Hence the inertia is primarily that of the pylon and rotor, with 

the wing contributing elastic restraint of the motion. The wing structural 

spring constants are , and ; these were evaluated by matching 

the predicted frequencies of^the wing modes to the values obtained experimentally 
* * * 




and C 

P 


are the structuxal damping constants for the wing modes. 


Vertical bending elevates the rotor trim thrust above the Inboard sections, 
and so gives a nose down pitch moment with effectiveness given by C^i 
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Dimensionally, the spring .'.nd n," constants are 

W*« k/$**tt* 

c* » e/« x w\si\ 

Hence the restive spring and damping rites vary with the rotor rotational 
speed} i.e. . ' wing frequency is really a fixed dimensional value (Hz), so 
the per-rev values vary with . 

Additional discussion and details of the wing equations of notion 
are given in reference 4. 

Wing Aerodynamics 

The wing aerodynamic forces exciting bending and torsion notion of 
the wing are: 
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where F and F are the vertical and chordwise aerodynamic forces on 
*w w t 

the wing section (lift and profile plus induced drag); is the aerodynamic 
noment about the local EA. The velocity seen by the section has perturbations 
due to the wing degrees of freedom, and due to aerodynamic gusts. Aerodynamic 
interference between the rotor and the wing is neglected . From the velocity 
perturbations, the perturbations of the section forces may be found, and hence 
the wing aerodynamic coefficients. The derivation of the wing aerodynamic 
coefficients follows the standard techniques of strip theory In aeroelasticity ; 
more details of the derivation are given in reference 4. We also include 
here the aerodynamic force due to the deflection of a control surface (flap 
or aileron) on the wing trailing-edge. The geometry is shown in figure 14, 

A constant chord (cp) trailing-edge flap, extending from y w * 
js considered. The flap deflection angle is , positive 


y F y w " y F 
downward • So ® 
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^ is a control variable, in addition to the rotor cyclic and collective 
pitch controls. The result for the wine aerodynamic forces is: 
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The aerodynamic coefficients arei 
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and C^ are the 

aircraft trim lift and drag (profile plus induced) 

coefficients} and 


and Cj^ their derivatives with respect to o<^ . 


The section monent characteristics are given by , the distance the wing 

AC is behind the EA, and c , the nose up moment Coefficient about the AC. 

in 

SIC 

The constant M 

- Cw 

nrra 

accounts for the difference ir. the normalization of the wing and rotor 
coefficients. The constants e R and f n are integrals of the wing mode shapes, 
accounting for the way the motion produces forces on the wingi 
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1 


For the flaperon coefficients, we usei 



, / , v C4 N , |-Cc /c«* v 

s w yl-t* 6 / 6^5 + **S c%* J 0 ^ 2/^ w } 


where 


c, * 


Z 


£/ 

C-w 


The first factor in these expressions is the two-dimensional thin airfoil 
theory result for the lift and moment due to control surface deflection; 
and the last two factors are corrections for the wing aspect ratio, thickness, 
and real flow effects on the flap effectiveness or- re^. (s'). 


tys Equation of Notion 

I 

The rotational speed decree of freedom ( *^ ) ic an important 
factor in the dynamics, especially with a windmilling rotor. Usually 
the equation of motion will involve the engine, drive train, interconnect 

shaft, and governor dynamics $ here we shall consider only two 1 Lilting cases. 


The first case is windmilling or autorotation operation of the 
rotor. The rotor is free to turn on the shaft, so no torque moments are 
transmitted from the rotor to the shaft, and no pylon roll motion transmitted 
to the rotor. Both effects are accomplished by using Cq * 0 as the 
equation of notion for Vg . There is no spring term on Vj , so the 
degree of freedom is really * the rotor speed perturbation. It should 
be noted that N*$ is defined with respect to the pylon, which has a roll 
angle o s so the rotor speed perturbation with respect to space is 
the sum 
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The second case considered here Is powered operation of the rotor. 
It Is assumed that the rotor hub rotational speed Is fixed, at n , with 
no perturbations. This case may be viewed as the limit of operation with 
a perfect governor on engine or rotor speed. The powered case is treated 
by dropping the degree of freedom and equation} i.e. the solution 

is just “ 0. 


Hence we add to the support equations of motion the equation 
for H'j j 




For the powered case this equation and the % degree of freedom are 
dropped from the system (a row and column eliminated from the matrices). 
For the windmilling case they are retained} note that the equation 
is first order, since there is no spring term. 

Reference 4 gives a further discussion of these two cases, 
windmilling and powered operation, and their effects on the pro pro tor 
dymaice. 


Support Equations of Motion 

We have obtained now the shaft motion and support equations of 
motion, which in natrix form arei 
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where the wing degrees of freedom ( x^ ) and the wing flap control 
( \ ) ere 
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v„ = C Sj ] 

and as defined above, the rotor hub forces and moments ( F ), shaft 
motion ( o< ), and aerodynamic gust ( 7^ ) arei 
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Ths na trices of the coefficients of the equations of motion follow. 
The aetrix Cj relating the rotor shaft notion to ths wing notion* has 
bssn given above. 
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EQUATIONS HCTION 

The conplete set of equations of motion describing the proprotor 
and cantilever wing system may now be obtained, by substituting for the 
shaft motion into the rotor forcer, and moments, and. then for the rotor forces 
Into the wing equations. , ''he result is a set of linear differential 
equations, of the forms 

A 2 x + A»x + Ao* — Bv + Q>& ^ 

where the degrees of freedom estate) vector ( "x ) and the input vector 



Recalling the equation® for the rotor equations of motion, the rotor hub 
forces and moments, the shaft notion, and the wing equations of motions 
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A|*g •+ AoX# 4 * A 2 oi + A| o* \ 0 a< — 8 v^ + 8»C<r^ 

F * C 2 %c -«* C, ** -*- c*** + Cj ♦ C, ©» + CoO< -f bft«j 
^ * ex«v 

'* t - 

«*» *w 4> Qo %w a 9> V«w -V •+» OS r 

the coefficient, matricer of the conrteto equation? of motion may be 
identified, or,: 



-141- 


i - i 


Treatment of Rotor Pltch/Torslon 

The equations of notion have been set up including the rotor pitch 
ami torsion degrees of freedom, © , and with © (the comnari^ed 

pitch angle) as the rotor control variable. One nay not wish to include 
these degrees of freedom in the system dynamics, but it is not possible to 
simply drop then at this stage. The pitch control and bend ing/ginbal 
feedback enters the system through the rigid pitch degree of freedom (p q )» 
so it is necessary to first operate on the columns of the equation matrices 
to account for these effects. Then the degrees of fi'eedom and equations 
(columns and rows of the matrices) nay be dropped as appropriate. We shall 
consider three options for the treatment of the rotor pitch/torsion motion. 


The first option is to include the pitch and torsion degrees of 
freedom in the system; then the equations are used as derived. 


The second option is the case of a rigid control system. It is 
the limit of infinite control system and blade torsion stiffness. Thus 
the rotor blade elastic torsion motion is zero, and the response of the 
rigid pitch motion reduces to 

(>o * e 1 *' - £ fcf i 
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Thus we operate on the columns of 

the 

A matrix as 

O 

follows : 

subtract Kp 

times 

the 

©r 

column 

from 

the 

i 



co> 




subtract K p 

times 

the 

e.c 

column 

from 

th 

i 

subtract K p 

times 

the 


column 

from 

the 

r i 







subtract Kp 

times 

the 

©ic 

column 

from 

the 

G 







subtract K p 

times 

the 

a* 

column 

from 

the 




(Vm 


column 

column 

column 

column 

column 


- 142 - 



and reconstruct the 
replace the 

replace the 

replace the 


control matrix B 


eT 

column of 

Gfrv* 

Su. 

column of 



column of 


as follows* 

B with minus the 

B with minus the 

)> with minus the 


<©> 

©o column of A 

o 

(o) 

column of A 

o 

.CO) 

column of A 

** o 


Then the rigid pitch degrees of freedom and equations of motion are 
dropped from the system. Note that the aoove transformation is only the 
result of infinite control system stiffness; 5t would be possible to retain 
the elastic torsion degrees of freedon, dropping only the rigid pitch p^. 


The third option is a quasistatic approximation for the effect 
of the blade torsion and pitch motion, ’.Je shall neglect the acceleration 
and velocity terms in the torsion/pitch equations. The torsion/pitch 
equations then become just a static substitution relation for Q in the 
other equations of motion. This treatment retains all the static couoling 
effects in the A matrix. The required transformation of the equations is 
accomplished as follows. First the A , B , and 3 matrices are partitioned, 

O Li 

to separate the ^ variables and equations from the rest. Assuming the 
© block is in the middle of x, the state equations take the forms 



Now the acceleration and velocity terms are dropped from the pitch equations; 
and ,e write x still for the state variable vector, but now with the 
pitch degrees of freedom dropped. Hence 


0 LA* A?]* + ^ j 


which may be substituted into the remaining equations, eliminating © from 
(the pitch acceleration and velocity terms in the remaining equations 
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are dropped). Thus the quasistatic torsion approximation giver the 
following equations of motion, in terms of the reduced state variable 
x (without the torsion/pitch degrees of freedom) t 
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Elastic ails (EA) 



figure ?. lotor hub and root geometry (undi.tr.rMd). showing gfarttl 
undersling (s i>A ). torque offset (xj, A ). feathering axis 
offset (r FA ), preeone < i*,) • droop (<»,) . -* e»eep ( *»A,> 
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Figure 9. Notation ami sign convention for gimbal motion, 
(a) in the nonrotating frame, and (b) in the 
rotating frame. 
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:odynamicsj notation and sign 
1 forces and velocities. 




Figure 11. Notation and Bits* convention* for rotor velocity and 
orientation (V and induced velocity (v), and 

aerodynamic gust velocity components (Uj,, v Q , « G ). 
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inure 14. Geometry of win," and wing flap aerodynamic 





